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Abstract 

We study the near-critical FK-Ising model. First, a determination of the correlation 
length defined via crossing probabilities is provided. Second, a striking phenomenon 
about the near-critical behavior of FK-Ising is highlighted, which is completely missing 
from the case of standard percolation: in any monotone coupling of FK configurations 
LOp {e.g., in the one introduced in |Gri95| ). as one raises p near pc., the new edges 
arrive in a fascinating self-organized way, so that the correlation length is not governed 
anymore by the number of pivotal edges at criticality. In particular, it is smaller than 
the heat-bath dynamical correlation length determined in the forthcoming [GP . 

We also include a discussion of near-critical and dynamical regimes for general 
random-cluster models. For the heat-bath dynamics in critical random-cluster models, 
we conjecture that there is a regime of values of cluster- weights q where there exist 
macroscopic pivotals yet there are no exceptional times. These are the first natural 
models that are expected to be noise sensitive but not dynamically sensitive. 



Contents 



1 Introduction [3] 



2 Proofs of the main results on the correlation length (Theorems 1.3, 1.2 



and[T]4]) [H 

2.1 Basic properties of random-cluster models 

2.2 Connectivity probabilities and the fermionic observable [13] 

2.3 Integr ability relations of the fermionic observable away from criticality . . . [TT] 

2.4 Proof of Theorem O l20l 



2.5 Proofs of Theorems 11.21 and 11.4 



Near-critical behavior: a fascinating self-organized near-criticality emerges 

3.1 Pivotal points govern the near-critical regime of percolation |26] 

3.2 Specific heat of the random-cluster model |29] 

3.3 The monotone increasing Markov process on cluster configurations |30] 

3.4 Existence of emerging clouds when q = 2 |34] 

3.5 Questions on the structure of emerging clouds |36] 

3.6 What about the influence of an edge? |37] 



4 Other values of q 

4.1 Some useful critical exponents |38l 

4.2 Near-critical behavior for q S (1,4] and self-organized monotone coupling . . 

4.3 The hyperscaling relation between correlation length and speciflc heat . . . [41] 

5 On the dynamical sensitivity of random-cluster models 1421 

5.1 Noise- versus dynamical sensitivity for g G [1,4] 22] 

5.2 Random-cluster models with q > 4 

References 



2 



1 Introduction 



Phase transition in the random cluster model. The random-cluster model with 
parameters p G [0, 1] and q > 1 [^is a probability measure on subgraphs of a finite graph 
G={V,E), defined for ah w C S by 

open edges ^-[^ closed edges^# clusters 

4>p,q{^) ■■= y , 

where Zp^q is the normalization constant such that (pp^g is a probability measure. The most 
classical example of the random-cluster model is bond percolation, which corresponds to 
the q = 1 case. Infinite volume measures can be constructed using limits of the above 
measures along exhaustions by finite subsets (with different boundary conditions: free, 
wired, etc). Random-cluster models exhibit a phase transition at some critical parameter 
Pc = Pc{q)- On Z"', this value does not depend on which infinite volume limit we are using, 
and, as in standard percolation, below this threshold, clusters are almost surely finite. 



while above this threshold, there exists (a.s.) a unique infinite cluster. See Subsection 2.1 
for details and references. 

The critical parameter is known to be equal to 1/2 for bond percolation on the square 
lattice. For the random-cluster model with parameter q = 2 (also called FK-Ising), Pc{2) = 
is known since Onsager |Ons44| (it is connected via the Edwards-Sokal coupling to the 
critical temperature of the Ising model). See also the recent |BDC10b| for an alternative 
proof of this fact. More recently, the general equality Pc{q) = "^^s proved for every 

g > 1 in [ BDClOa] . 

Mathematicians and physicists are interested in the properties of the critical phase 
{p,q) = {Pc{q)^q) itself. Very little is known for general values of q, and only q = 1 and 
q = 2 are understood in a satisfactory fashion. For these two cases, the phase transition 
is known to be of second order [i.e. continuous): at pc, almost surely there are only finite 
clusters. For q = 2, Smirnov proved the conformal invariance of certain macroscopic 
observables |SmilO| ISmi| . which, combined with the Schramm-Lowner Evolution |SchOO| . 
can be used to construct continuum scaling limits that retain the macroscopic cluster 
structure. For q = 1, site percolation on the triangular lattice (a related model) has been 
proved to be conformally invariant as well |Smi01| . 

Beyond the understanding of the critical phase, the principal goal of statistical physics is 
to study the phase transition itself, and in particular the behavior of macroscopic properties 
(for instance, the density of the infinite-cluster for p > Pc{q))- It is possible to relate the 
critical regime to these thermodynamical properties via the study of the so-called near- 
critical regime. This regime was investigated in |Kes87| in the case of percolation. Many 
works followed afterward, culminating in a rather good understanding of dynamical and 
near-critical phenomena in standard percolation |GPS10al IGPSlOb"! IGPSj . The goal of this 
article is to discuss the near-critical regime in the random-cluster case, and more precisely 
in the FK-Ising case. 



Near-critical regime and correlation length. The near-critical regime is the study 
of the random-cluster model of edge-parameter p in the box of size L when (p, L) goes to 
(pc) oo). Note that, on the one hand, if p goes to Pc very quickly the configuration in the 

^In general one could take g > 0, but we will assume q>l here, in order for the FKG inequality to 
hold, see |Gri06| . 
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box of size L will look critical. On the other hand, if p goes to Pc (from above) too slowly, 
the random-cluster model will look supercritical. The typical scale L = L{p) separating 
these two regimes is called the correlation length (or characteristic length). In rough 
terms, if p is slightly above Pc(2) = ^2/(1 + \/2), the correlation length L{p) is the scale 
below which things still look somewhat critical and above which the infinite cluster starts 
being visible. In the subcritical regime, it corresponds to the scale above which the fact 
that p is subcritical becomes apparent. 



Definition of correlation length in the case of percolation (g = 1). The critical 
regime is often characterized by the fact that crossing probabilities remain strictly between 
and 1. Formally, consider rectangles R of the form [0,n] x [0, m] for n,m > 0, and 
translations of them. We denote by Cy{R) the event that there exists a vertical crossing 
in R, a path from the bottom side [0, n] x {0} to the top side [0, n] x {m} that consists 
only of open edges. The celebrated Russo-Seymour- Welsh theorem shows that in the case 
of critical percolation, crossing probabilities of rectangles of bounded aspect ratio remain 
bounded away from and 1. A natural way of describing the picture as being critical is to 
check that crossing probabilities are neither near nor near 1. Mathematically, it is thus 
natural to define the correlation length for every p < pc = 1/2 and e > as 

Ls{p) := mf{n>0 : Fp{C^{[0,nf)) <e} , 

and, when p > pc = 1/2, as Li,{p) := L^{1 — p), where 1 — p is the dual edge-weight. The 
dependence on e is not relevant, since it can be proved f |Kes871 iNoM] ) that for any e > 0, 

^1/4 (P) - Le{p) , 

where x means that there exist constants < As, < oo such that 

A,Ly^{p) < Le{p) < BsLy^ip). 

The correlation length was shown to behave like \p — pc{l)\~'^^^~^°^^^ in the case of 
percolation |SW01j . 



Correlation length for FK-Ising percolation {q = 2). The first result of this paper 
is a determination of the behavior of L{p) when p goes to Pc for q = 2. Before stating 
the main result, let us give a proper definition of correlation length in this setting. Since 
the Russo-Seymour- Welsh theorem has been generalized to the FK-Ising case in [DCHNIO] 



(see Theorem 2.1 in the present text), it is natural to characterize the critical regime once 
again by the fact that crossing probabilities remain strictly between and 1. An important 
difference from the q = 1 case is that one has to take into account the effect of boundary 



conditions (see Subsection 2.1 for precise definitions): 

Definition 1.1 (Correlation length). Fix q = 2 and p > 0. For any n > 1, let Rn be the 

rectangle [0, n] x [0,/3n]. 

If p < Pc{2), for every e > and boundary condition ^, define 

L%,{p) := mi {n > : cj^l^^^JCiRn)) < e} , 

where (t>\ip2 denotes the random- cluster measure on R with parameters {p, 2) and boundary 
condition S^. If p > Pc{2), define similarly 

LiM ■■= inf {n > : ^ 2,r„ [C,{Rn)) >l-e}. 
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Our main result on the correlation length can be stated as follows. 




Theorem 1.2. Fix q = 2. For every e,p > 0, there is a constant c = c{€, p) such that 

\P - Pel ^' 

for all p 7^ pc whatever the choice of the boundary condition ^ is. 

Note that the left-hand side of the previous theorem has the following reformulation, 
which we state as a theorem of its own (this result is interesting on its own since it provides 
estimates on crossing probabilities which are uniform in boundary conditions away from 
the critical point): 

Theorem 1.3 (RSW-type crossing bounds). For A > and p > 1, there exist two constants 
< c_ < c+ < 1 such that for any rectangle R with side lengths n and m € [^n, pn], any 

p G [Pc — ^iPc -\- ^ (ind any boundary condition ^, one has 

C-<4,p,2(C.(^))<C+. 

The main ingredient of the proof of the latter theorem (and the most interesting one) is 
Smirnov's fermionic observable. This observable is defined in Dobrushin domains (with a 
free and a wired boundary arc) , and is a key ingredient in the proof of conformal invariance 
at criticality. Nevertheless, its importance goes far beyond that proof, in particular because 
it can be related to connectivity properties of the FK-Ising model. We study its properties 
away from the critical point (developing further the methods of |BDC10b] ). and estimate 
its behavior near the free arc of Dobrushin domains. It implies estimates on the probability 
for sites of the free arc to be connected to the wired arc. This, as in [DCHNIO] . allows us 
to perform a second-moment estimate on the number of connections between sites of the 
free arc and the wired arc, therefore implying crossing probabilities in Dobrushin domains. 
All that remains is to get rid of the Dobrushin boundary conditions (which is not as simple 
as one might hope; in particular, harder than in [DCHNlO] ). in order to obtain crossing 
probabilities with free boundary conditions. 

Using conformal invariance techniques, Chelkak, Hongler and Izyurov has recently 
proved that 

(f>p^,q=2{0^d[-n,nf) ~ Cn-'/\ (1.1) 

together with the appropriate (conformally invariant) version for general domains |CHI12) . 
Together with Theorem |1. 2 [ this implies: 



Theorem 1.4. Assuming (1.1), there exists a constant c > such that if p > Pc(2), 

\ 1/8 



</'p,2(0 -H- OO) > C 



\P-Pc 



nogi/\p-pc\ J 

The critical 1-arm exponent ( |1.1[ ) and the off-critical result 

{cTo)^ = 0p,2(O oo) X 1/3 - /3,|i/8 (as p \ P,) 
go back to Onsager |Ons44| . These two results more-or-less imply the correlation length 

\P-Pc\ 
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with more precise and direct arguments presented in |FF69| and |Kad69| (see also [BDClOb] ). 
However, their notion of correlation length is different and more restricted than the one 
in our Theorem 1.2, and its equivalence with our notion is not known. Furthermore, 



our 



proof of Theorem 1.4 is also of some value, since the result of |CHI12j and the techniques 
in this paper extend to isoradial graphs (with additional work) while Onsager's technology 
is restricted to the square lattice. 



The random-cluster model through its phase transition. The previous way to 
look at the near-critical regime may seem slightly artificial. It is more natural to study the 
random-cluster model through its phase transition by constructing a monotone coupling of 
random-cluster models with fixed cluster- weight q > 1. Then, properties of the monotone 
coupling (which can be thought of as a dynamics following the evolution of p between 
and 1) near pc will describe the near-critical regime. 

In the case of standard bond percolation {q = 1), such a monotone coupling simply 
consists of i.i.d. Uniform[0, 1] labels on the edges, and a percolation configuration LVp of 
density p is the set of bonds with labels at most p. It is straightforward to interpret 
this coupling as an asymmetric dynamical percolation: starting from critical percolation 
at time zero, as time goes on, whenever the clock of a bond rings, we open that bond; 
we can also run time backwards and close the bonds that ring. Now, the question is: in 
this monotone coupling, how fast does the system enter the supercritical and subcritical 
regimes as p changes near 

This near-critical window in percolation was studied by Kesten in |Kes87| . then by 
|B(:;KSni| INoinSl INWn9| l(;tPSinbl EEs], it turns out that its size is governed by the 
expected number of macroscopically pivotal edges at criticality, i.e., edges having four 
alternating (between dual and primal) open paths starting there and going to a macroscopic 
distance. Let us describe briefly how this mechanism works, since one of our main goals 
in this article is to point out that, in the case of FK-percolation (g > 1), a striking new 
phenomenon appears. 

Let 04 (n) be the probability at criticality that an edge has four alternating paths going 
to distance n. Getting from ujp^ to Up^+Ap in the box of size L, the system is moving out of 
stationarity, and roughly L'^Ap edges are switched from closed to open (we are assuming 
Ap > 0). The expected number of opened edges that were closed macroscopic pivotals 
in the initial configuration (preventing macroscopic open paths) is about L^a4(L)Ap. If 
L?'ai{L)Ap ^ 1, then many of these initial macroscopic pivotals have become open with 
good probability. It is not very hard to show that this means that the window of size L 
has become well-connected, hence we have left the near-critical regime. 

What can we say if L?'a/i{L)Ap <^ 1? The number of initial macroscopic pivotals 
that have switched is small, but maybe many new pivotals have appeared during the 
dynamics, which could have switched then, establishing macroscopic open connections. To 
formulate the same issue from a different point of view, if the dynamics, instead of switching 
always from closed to open, was symmetric dynamical percolation, then the system would 
be critical all the time, and, using Fubini and the linearity of expectation, the expected 
number of macroscopic pivotal switches (i.e., flips of edges that are macroscopically pivotal 
at the moment of the flip) in time Ap would be L^a4(L)Ap. If this expectation is small, 
then the probability of having any macroscopic pivotal switches is also small, hence the 
system indeed has not changed macroscopically. However, the asymmetric near-critical 
dynamics is slowly moving out of criticality, which could have an effect on the number of 
pivotals, speeding up changes. 
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Nevertheless, Kesten proved the following near-critical stability result: as long as 
L'^a4{L)Ap = 0(1), there are not many more pivotal points in oOp than at criticality, hence, 
despite the monotonicity of the dynamics, changes do not speed up significantly compared 
to symmetric dynamical percolation, and hence the macroscopic geometry starts changing 
significantly only at L^a4(L)Ap x 1. The proof in |Kes87| employs Russo's influence 
formula and differential inequalities, discussed more in Subsection |3.1| There is a related 
but more geometric approach in |GPS| . which provides some additional insight into the 
mechanism governing the near-critical regime that will be crucial in understanding the case 
of FK-percolation. 

Namely, (GPSj proves the following dynamical stability result about symmetric dy- 
namical percolation (needed to establish that dynamical percolation has a scaling limit that 
is a Markov process describing the changes of the macroscopic connectivity structure): as 
long as L^a4(L)At = 0(1), in order to describe the macroscopic structure of WAti it is 
enough to know the macroscopic structure of ujq and to follow the flips experienced by all 
initial macroscopic pivotals. In other words, there are no cascades of information from 
small to much larger scales, i.e., edges initially pivotal only on a "mesoscopic" scale are 
unlikely to have a macroscopic impact in the dynamics within the given time frame. This is 
proved using induction, with a careful summation over all possible ways in which "smaller" 
pivotals can make a big difference. For this summation argument, it is of course essential 
(and this will be important later) that edges are switching independently from each other. 
Now, the same argument applies to the asymmetric dynamics, and gives the near-critical 
stability that we stated above: in order to change the macroscopic connectivity structure, 
initial macroscopic pivotals need to be flipped. 

Summarizing, the scale at which the critical regime becomes the supercritical regime 
is given by L'^ai{L)Ap x 1. The same reasoning can be applied to the subcritical regime. 
In particular, the correlation length is given (up to constants) by Kesten's relation 

{Ls{p))^a4{L,{p))\p-Pc\ ^ 1, (1.2) 

where the constant factors of x depend only on the value of e. 

Going back to dynamical percolation for a second, it is not very hard to show that 
L^Q;4(L)At > c > implies that w^t is already somewhat different from loq, while the 
much harder [GPSlOa] says that L'^a4^{L)At ^> 1 implies that the macroscopic structure 
of has completely forgotten the initial configuration lvq. But even without the latter 
result, we see that, for any time interval At, if L = L{At) is given by L'^a^{L)At x 1, 
then, below this scale L, the macroscopic structure is essentially untouched, while, at larger 
scales, changes are already visible. So, this L may be called the dynamical correlation 
length. 

The main principle we shall extract from this discussion can be stated as follows: 

Phenomenon 1.5. In percolation (q = 1), due to near-critical stability, the near-critical 
behavior is governed by the number of pivotal points at criticality. The near-critical corre- 
lation length coincides with the dynamical correlation length. 

To our knowledge, it has been widely believed in the community that basically the same 
mechanism should hold in the case of random-cluster models. Namely, once we understand 
the geometry of the set of pivotal points at criticality, we may readily deduce information 
on the near-critical behavior. However, this is not the case. 

Let us consider the case of the FK-Ising. It is shown in |DCG) that the critical FK- 
Ising probability af^{n) for a site to be pivotal behaves like 7),-35/24+o(i) ^i^^-j^ ^ goes to 
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infinity. If pivotal points were governing the near-critical regime, the correlation length 
should satisfy 

{Lf^{p)fa^{Ll'^{p))\p-p,\ X 1, (1.3) 

which would give 

^P(P) = |p-Pc(2)|-^+°« » \p-p,{2)\-\ (1-4) 



contradicting Theorem 1.2 



What goes wrong with the percolation arguments? First, there is a basic phenomenon 
in the FK(p, q) models for q >2 that is very relevant to the above discussion: the difference 
between the average densities of edges for p = Pc{q) + and p = Pc{q) is not proportional 
to Ap, but larger than that, with an exponent given by the so-called specific heat of the 



model. (We will discuss this in more detail in Subsection 3.2 ) A first guess could be that 



the discrepancy in (1.4) is a result of the fact that Ap is not the density of the new edges 
arriving, and this should have been taken into account in the computation using the pivotal 
exponent. However, this is only partially right: the specific heat exponent itself is not large 
enough to account for this discrepancy (in fact, for g = 2 it equals — there is only a 
logarithmic blow-up). In fact, a self-organizational mechanism kicks in. In standard 
percolation, on the way from ojp^ to Wp^+Ap in the asymmetric dynamics, new edges arrive 
in a "Poissonian" way. This is no longer the case with FK-Ising: the arriving edges tend 
to prefer "strategic" locations, which are pivotal at large scales, thereby speeding up the 
dynamics compared to what could be guessed from the number of pivotals at criticality. 
In other words, near the arriving edges depend in a very sensitive way on the current 
configuration. This subtle balance between the current configuration and the conditional 
law of the arriving edges is representative of a self-organized mechanism. 

How do we know about this self-organization? Similarly to dynamical percolation, 
there is a natural dynamics with the random-cluster model FK(p, q) as stationary distri- 



bution, called the heat-bath dynamics or Sweeny algorithm (see Definition 3.3): edges have 
independent exponential clocks, and when the clock of e = (x, y) rings, the state of e is 
updated according to the FK(p, q) measure conditioned on the rest of the configuration w, 
of which the only relevant information is whether x and y are connected in w \ {e}. Now, 
in |GPj . the analogue of the above-mentioned dynamical stability result of |GPS| is proved 
for the critical FK-Ising model FK(pc(2), 2). Thus, if there was any monotone coupling of 
the near-critical FK(p, 2) models in which new edges arrived in a Poissonian way similar to 
the heat bath dynamics, with clock rates and resampling probabilities bounded away from 
0, then the same argument would apply, and near-critical stability would hold, proving 



(1.3). Since this is not the case, the Poissonian picture can be ruled out: 



Phenomenon 1.6. The correlation length in FK-Ising is much smaller than what the 
intuition coming from standard percolation (q = 1) would predict. In other words, as one 
raises the parameter p, the supercritical regime appears "faster" than what would he dictated 
simply by the number of pivotal edges at criticality: new edges arrive in a very non-uniform 
manner, and a striking self-organized near-criticality appears. 

The proof that this self-organization mechanism must play an important role in any 
monotone coupling depends on the forthcoming dynamical FK-paper |GP| . but we never- 
theless wanted to draw attention to this very surprising near-critical corollary already in 
this paper. A particular realization of this phenomenon that we can offer here is the follow- 



ing. In Subsection 3.3, we will introduce Grimmett's monotone coupling of the FK(p, ( 
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configurations as p varies from to 1 and g > 1 is fixed. In Subsection 3.4, we will 
present a concrete way in which self-organization works in this coupling. However, most 
of the self-organization scheme remains to be understood; some questions are presented 
in Subsection |3.5| Furthermore, it is interesting to note that there is a hyperscaling rela- 
tion between specific heat and correlation length, hence the specific heat determines not 
only the number of new edges arriving, but also the correlation length (in some way that 
appears to be independent of the self-organization scheme). We will briefly discuss this 
hyperscaling in Subsection |4.3[ 



In Subsection 3.6 we present another point of view that might explain the discrepancy 
(1.4): for q > 1, Russo's formula used in Kesten's proof has to be modified. Namely, the 
influence of an edge on the event that a box of size n is crossed does not coincide with 
the probability for that edge to be pivotal, as it was the case for q = 1- Let us define the 
critical exponent by assuming that the above influence of an edge behaves like n~^^'^^ 
at criticality. Kesten's scaling relation for q = 1 was (2 — (,4{q))i'{q) = 1, where ^4(9) and 
i'{q) are the critical exponents of the pivotal event and the correlation length, respectively. 



coming from (1.2). This will remain valid for q > 1 only if the critical exponent S,4{q) is 
replaced by the exponent L{q) governing the behavior of the influence. This subtlety and 
the fact that ^4((?) 7^ i{q) seem to be new. 

On the dynamical and near-critical behavior for other values of q. Finally, we 
investigate what happens for other values of q. Since the mathematical understanding of 
critical FK{q) models is very limited when q {0,1,2}, the study relies on predictions 
from physics. We will mostly focus on the case q G [1,4], where the FKG inequality holds 
and the phase transition is conjectured to be continuous {i.e., there is a unique infinite- 
volume measure at criticality, having no infinite cluster). It is therefore natural to consider 
the near-critical regime. In this case yet again, pivotal points do not seem to control the 
behavior of the near-critical regime. Critical exponents indeed violate Kesten's relations. 
See Section lU 

We will also discuss briefly noise-sensitivity and dynamical sensitivity of random-cluster 
models with q £ [1,4]. Note that, as it is already apparent from the above discussion, 
the study of the near-critical regime of percolation (especially in [GPSlObl IGPS] ) was 
conducted in parallel to the study of dynamical percolation |SS10| IGPSlOa] . For random- 
cluster models with q E (1,4], one may indeed study the influence of pivotal edges on 
the existence of exceptional times in the heat-bath dynamics: times for which an infinite 
cluster exists. As for the near-critical regime, the situation seems much more complicated 
than the percolation one, although for different reasons. (For dynamical sensitivity, instead 
of the true near-critical model, a "fake, Poissonian" near-critical model appears to be more 
relevant.) In particular, the existence of pivotal edges is not equivalent (conjecturally) to 
the dynamical sensitivity of the model (it is expected to be equivalent to noise sensitivity. 



though). We refer to Subsection 5.1 for further details on these phenomena. 

Random-cluster models with q > 4 do not have the same behavior as those with g < 4. 
Indeed, the phase transition is conjectured to be of first order, and the near-critical regime 
exists only in a restricted sense: there is a non-trivial critical window in which a finite 
system becomes supercritical from subcritical, but the correlation length remains bounded 
for all p. In addition, the models are not expected to be noise-sensitive or dynamically 



sensitive. Subsection 5.2 is devoted to their study. 
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Organization of the paper. In Section [2] we begin with an exposition of basic prop- 
erties of random-cluster models, and define Smirnov's so-called fermionic observable. We 



use it away from criticality to prove Theorem |1.3[ then we deduce Theorems 1.2 and 1.4 



In Section [3j we study the self-organized phenomenon in detail. We start by explaining 
why pivotal points are crucial in the understanding of the near-critical regime of percola- 
tion. First, we discuss the role of specific heat, then present Grimmett's monotone coupling, 
under which the evolution of the random-cluster model through its phase transition can 
be monitored. Finally we explain how the self-organized phenomenon acts concretely. 

Section |4] contains a discussion of other values of while Section [5] mentions the 
interesting case of dynamical random-cluster models. 



2 Proofs of the main results on the correlation length (The- 
orems 



1.3, 1.2 and 1.4) 



In this section, a point will be identified with its complex coordinate. 
2.1 Basic properties of random-cluster models 

The random-cluster measure can be defined on any graph. However, we restrict ourselves 
to the standard square lattice Z^. With a tiny abuse of notation, we will use V{1?) or just 
1? for the set of sites, and E{7?) for the set of bonds. In this paper, G will always denote 
a connected subgraph of Z^, i.e., a subset of vertices together with all the bonds between 
them. We denote by dG the (inner) boundary of G, i.e., the set of sites of G linked by a 
bond to a site of I? \ G. 

A configuration w on G is a random subgraph of G, having the same sites and a subset 
of its bonds. We will call the bonds belonging to oj open, the others closed. Two sites a 
and h are said to be connected (denoted by a -f-T- 6), if there is an open path — a path 
composed of open bonds only — connecting them. The (maximal) connected components 
will be called clusters. More generally, we extend this definition and notation to sets in a 
straightforward way. 

A boundary condition ^ is a partition of dG. We denote by w U ^ the graph obtained 
from the configuration lo by identifying (or wiring) the vertices in ^ that belong to the same 
class of A boundary condition encodes the way in which sites are connected outside of 
G. Alternatively, one can see it as a collection of abstract bonds connecting the vertices in 
each of the classes to each other. We still denote by cj U ^ the graph obtained by adding 
the new bonds in ^ to the configuration oj, since this will not lead to confusion. Let o{ijj) 
(resp. c{(jj)) denote the number of open (resp. closed) bonds of cj and /c(a;,.^) the number 
of connected components of a; U ^. The probability measure (/)^ p g random-cluster 
model on a finite subgraph G with parameters p G [0, 1] and q G (0, oo) and boundary 
conditions ^ is defined by 

(M) := ^ ^ ^ , (2.1) 



G,p,g 

for any subgraph oj of G, where ^cpg ^ normalizing constant known as the partition 
function. When there is no possible confusion, we will drop the reference to parameters in 
the notation. 
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The domain Markov property. One can encode, using an appropriate boundary con- 
dition ^, the influence of the configuration outside a sub-graph on the measure within it. 
Consider a graph G = iy,E) and a random-cluster measure (Pcpg it- For F C E, 
consider G' with F as the set of edges and the endpoints of it as the set of sites. Then, the 
restriction to G' of (j)Q ^ ^ conditioned to match some configuration ui outside G' is exactly 
(f>Q p , where ^ describes the connections inherited from a; U (two sites are wired if they 
are connected by a path in u Dip outside G' — see Lemma 4.13 in |Gri06] ). This property 
is the direct analog of the DLR conditions for spin systems. 

Comparison of boundary conditions when q > 1. An event is called increasing if 
it is preserved by addition of open edges. When q > 1, the model satisfies the FKG- 
inequality, or is positively associated (see Lemma 4.14 in |Gri06| ). which has the following 
consequence: for any boundary conditions ip < ^ (meaning that ip is finer than ^, or in 
other words, that there are fewer connections in ip than in ,^), we have 

ApM) ^ AvM) (2.2) 

for any increasing event A. This last property, combined with the domain Markov property, 
provides a powerful tool to study the decorrelation between events. 

Examples of boundary conditions: free, wired, Dobrushin. Three boundary con- 
ditions play a special role in the study of random-cluster models: 

• The wired boundary conditions, denoted by (pQpq, is specified by the fact that all 
the vertices on the boundary are pairwise connected. 

• The free boundary conditions, denoted by (pQ ^ is specified by the absence of wirings 
between boundary sites. 

These boundary conditions are extremal for stochastic ordering, since any boundary con- 
dition is smaller (resp. greater) than the wired (resp. free) boundary conditions. 

• The Dobrushin boundary conditions: Assume now that dG is a self-avoiding polygon 
in L, let a and b be two sites of dG. The triple {G, a, b) is called a Dobrushin domain. 
Orienting its boundary counterclockwise defines two oriented boundary arcs dab and 
dfja] the Dobrushin boundary conditions are defined to be free on dab (there are no 
wirings between boundary sites) and wired on dba (all the boundary sites are pairwise 
connected) . These arcs are referred to as the free arc and the wired arc, respectively. 
The measure associated to these boundary conditions will be denoted by fpcpq o^' 
simply (/»^^ 

Infinite-volume measures and the definition of the critical point. The domain 
Markov property and comparison between boundary conditions allow us to define infinite- 
volume measures. Indeed, one can consider a sequence of measures on boxes of increasing 
sizes with free boundary conditions. This sequence is increasing in the sense of stochas- 
tic domination, which implies that it converges weakly to a limiting measure, called the 
random-cluster measure on with free boundary condition (and denoted by (pp^g). This 
classic construction can be performed with many other sequences of measures, defining 
several a priori different infinite-volume measures on Z^. For instance, one can define the 
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random-cluster measure ^ with wired boundary condition, by considering the decreasing 
sequence of random-cluster measures on finite boxes with wired boundary condition. 

On Z"^, for a given q > 1, it is known that uniqueness of the infinite- volume measure can 
fail only for p in a countable set Dq, see Theorem 4.60 of |Gri06j . Since all limit measures 
are sandwiched between cp^ ^ and (j)^ ^ w.r.t. stochastic domination, the countability of Vg 
implies that there exists a critical point pc such that for any infinite-volume measure with 
p < Pc (resp. p > Pc), there is almost surely no infinite component of connected sites (resp. 
at least one infinite component). 

Planar duality. In two dimensions, a random-cluster measure on a subgraph G of Z'^ 
with free boundary conditions can be associated with a dual measure in a natural way. 
First define the dual lattice (Z^)*, obtained by putting a vertex at the center of each face of 
7?, and by putting edges between nearest neighbors. The dual graph G* of a finite graph 
G is given by the sites of (Z'^)* associated with the faces adjacent to an edge of G. The 
edges of G* are the edges of (Z^)* that connect two of its sites - note that any edge of G* 
corresponds to an edge of G. 

A dual model can be constructed on the dual graph as follows: for a percolation 
configuration cj, each edge of G* is dual-open (or simply open), resp. dual-closed, if the 
corresponding edge of G is closed, resp. open. If the primal model is a random-cluster 
model with parameters {p,q), then it follows from Euler's formula (relating the number 
of vertices, edges, faces, and components of a planar graph) that the dual model is again 
a random-cluster model, with parameters {p*,q*) - in general, one must be careful about 



the boundary conditions. For instance, on a graph G, the random-cluster measure (j)Q. 
is dual to the measure (pQ* p* , where {p* , q* ) satisfies 

pp* 



p,<i 



(1 — p){l — p* 



q and q 



Similarly, the dual of Dobrushin boundary conditions are Dobrushin boundary conditions 
themselves. 

The critical point Pc{q) of the model is the self-dual point Psd(Q) for which p = p* (this 
has been recently proved in |BDC10a] ). whose value can be derived: 



1 + 



In the following, we need to consider connections in the dual model. Two sites x and 
y of G* are said to be dual- connected if there exists a connected path of open dual-edges 
between them. Similarly to the primal model, we define dual clusters as maximal connected 
components for dual-connectivity. 



FK-Ising model: crossing probabilities at criticality. For the value g = 2 of the 
parameter, the random-cluster model is related to the Ising model. In this case, the 
random-cluster model is now well-understood. The uniqueness of the infinite volume limit 
for all p is known since Onsager; see |Wer09b] for a short and elegant proof (or (DCSIH 
Proposition 3.10] for a version of Werner's proof in English). The value pc = Psd is 
implied by the computation by Kaufman and Onsager of the partition function of the 
Ising model, and an alternative proof has been proposed recently by Beffara and Duminil- 
Copin iJBDClObj. Moreover, in |SmilO| . Smirnov proved conformal invariance of this model 
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at the self-dual point Ps^- As mentioned in the Introduction, criticality can sometimes be 
characterized by the fact that crossing probabilities are bounded uniformly away from 
and 1. This fact was proved in jDCHNlO] (our result is an extension of this one away from 
criticality) : 

Theorem 2.1 (RSW-type crossing bounds, |DCHN10] ). Let < pi < p2- There exist two 
constants < c_ < c+ < 1 (depending only on p\ and P2) such that for any rectangle R 
with side lengths n and m G [/0l?^, P2'n\ ('i-e. with aspect ratio hounded away from and 00 
hy pi and P2), one has 

for any boundary conditions ^. 

In the rest of this section, we consider only random-cluster models on the 
two-dimensional square lattice with parameter q = 2, hence we drop the dependency 
on q in the notation. In this case, the model is called FK-Ising model. 



2.2 Connectivity probabilities and the fermionic observable 

The random-cluster with cluster-weight (7 = 2 is a model with long-range dependence. In 
particular, boundary conditions play a crucial role in connectivity probabilities. While 
general percolation arguments are sometimes sufficient to estimate connectivity probabil- 
ities in the bulk |BDC10a] . there are very few possibilities to control probabilities in the 
presence of boundary conditions. We thus need a new argument to control these crossing 
probabilities. 

When q = 2, Smirnov's fermionic observable provides us with a powerful tool to study 
such probabilities. In the next paragraph, we introduce the loop representation of the 
random-cluster model and we define Smirnov's observable. In the next one, we remind 
several properties of this observable at and away from the critical point. We list them 
without proof, since they are already presented in various places. 



The medial lattice and the loop representation. Let G be a finite subgraph of 
together with Dobrushin boundary condition given by the boundary points a and b. Let 
G* be the dual graph, with the natural definition that respects the boundary condition, see 
the left side of Figure 2.1 Declare black the sites of G and white the sites of G* . Replace 



every site with a colored diamond, as in the right side of Figure 2.1 The medial graph 



Go = (^j-E'o) is defined as follows: is the set of diamond sides which belong to both 
a black and a white diamond; is the set of all the endpoints of the edges in We 
obtain a subgraph of a rotated (and rescaled) version of the usual square lattice. We give 
Go an additional structure as an oriented graph by orienting its edges clockwise around 
white faces. 

The random-cluster measure on {G, o, b) with Dobrushin boundary conditions has a 
rather convenient representation in this setting. Consider a configuration co. It defines 
clusters in G and dual clusters in G*. Through every vertex of the medial graph passes 
either an open bond of G or a dual open bond of G*, hence there is a unique way to draw 
Eulerian {i.e., using every edge exactly once) loops on the medial lattice — interfaces, 
separating clusters from dual clusters. Namely, a loop arriving at a vertex of the medial 
lattice always makes a 7r/2 turn so as not to cross the open or dual open bond through 
this vertex, see Figure [2]T] Besides loops, the configuration contains a single curve joining 
the vertices adjacent to a and b, which are the only vertices in with three adjacent 
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wired arc 



9 ^--9 

; 

9 I o^-9---9--^ 

9--^--9--^---9 



a 



9 94-9---9--^ 
o---o4-9-^-9---o 

I 9 9 9 [ free arc 

6 - - -o - - -o - - -o 




Figure 2.1: Left: A graph G with Dobrushin boundary conditions, and its dual G* . The 
black (resp. white) sites are the sites of G (resp. G*). The open bonds of G (resp. G*) are 
represented by solid (resp. dashed) black bonds. Right: Construction of the medial lattice 
and the loop representation: the loops are interfaces between primal and dual clusters. 



edges. This curve is called the exploration path and is denoted by 7. It corresponds to the 
interface between the cluster connected to the wired arc and the dual cluster connected to 
the free arc. 

The first and last edge of 7 are denoted by Ca and eb, respectively. More generally, Cc 
denotes the medial edge pointing north-east and bordering the diamond associated to c. 

This construction gives a bijection between random-cluster configurations on G and 
Eulerian loop configurations on G^. The probability measure can be nicely rewritten 
(using Euler's formula) in terms of the loop picture: 



"(p)^ open bonds y/2' 

Z{p,G) 



# loops 



where x{p) :- 



P 



[1-pW2 



and Z{p, G) is a normalizing constant. Notice that p = pc ii and only if x{p) = 1. This 
bijection is called the loop representation of the random-cluster model. The orientation of 
the medial graph gives a natural orientation to the interfaces in the loop representation. 



The fermionic observable. Fix a Dobrushin domain (G, a,6). Following |SmilOj . we 
now define an observable F on the edges of its medial graph, i.e. a function F : £"<> — t- C. 
Roughly speaking, i*" is a modification of the probability that the exploration path passes 
through a given edge. 

First, the winding Wr{z,z') of a curve P between two edges z and z' of the medial 
graph is the total rotation (in radians and oriented counter-clockwise) that the curve makes 
from the mid-point of edge z to that of edge z' . We define the observable F = Fp for any 
edge e S £"0 as 

He) ■■= cp;ia (e'^'^^^^'^'K^,) , (2.3) 
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where 7 is the exploration interface from a to b. 



Relation with connectivity probabilities. As was mentioned earher, the fermionic 
observable is related to connectivity properties of the model via the following fact: 

Lemma 2.2 (Equation (14) in |SmilOj . Lemma 2.2 in |BDC10b] ). Let u £ G be a site 
next to the free arc, and e be a side of the black diamond associated to u which borders a 
white diamond of the boundary. Then, 

|F(e)| = 4'p'2g('^ ^ wired arc). (2.4) 

Integrability relations of the fermionic observable. A lot of information has been 
gathered on the fermionic observable during the last few years. In particular, it satisfies 
local relations that allow to determine its scaling limit. 

Proposition 2.3 (Lemma 2.3 of |BDC10bj ). Consider a medial vertex v in G^\ dG^ . We 

index the two edges pointing towards v by A and C, and the two others by B and D, such 
that the alphabetical order is clockwise oriented. Then, 

F{A)-F{C) = i[F{B) - F{D)]. (2.5) 

where 



+ X 



(2.6) 



Moreover, the complex argument modulo vr of at any edge e is determined by the 
direction of e. More precisely, if e points in the same direction as Ch, then F is real. 
Similarly, F belongs to e~*'^/^M (resp. iM, e^'^/^M) if e makes an angle of tt/2 (resp. vr, 
37r/2) with the edge €},, see Fig. |2.2[ Knowing the complex argument modulo vr, together 



with (2.5), allow one to express the value of the observable at one edge e in terms of the 
values at two edges incident to one of the endpoints of e. This important fact was used 
extensively in |SmilO| and in any following work since. For instance, it implies that F is 



determined by the relations (|2.5[) and the fact that F{eb) = 1 and F{ea) = e*^r(eo,e6)/2^ 



where T is any path from to Ch staying in G^. In addition, these relations have a very 
special form. In particular, it has been proved for x = 1 in. |SmilO| and then extended 
in |BDC10b| (Lemma 4.4) to general values of x that F is massive harmonic inside the 
domain: 

Proposition 2.4. Let p € (0, 1) and X with four neighbors in G\ dG, we have 

ApF{ex) = 0, (2.7) 

where the operator Ap is defined by 



Apgiex) := V g{ey) - g{ex). (2.8) 




Observe that a{p) = if and only if p = pc. In this case, the observable is discrete 
harmonic inside the domain. As mentioned before, this is one of the main ingredients 
of Smirnov's proof of conformal invariance: when properly rescaled, the observable con- 
verges to a harmonic map. Boundary conditions for F correspond to discretizations of the 
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Figure 2.2: Left: An edge inside the domain: it has four edges oriented the same way at 
distance two. Right: An edge on the free arc with the associated indexation. 

Riemann-Hilbert problem. These boundary conditions are quite comphcated to study at 
a discrete level, and Smirnov used a discrete primitive H of (the imaginary part of) F'^ to 
handle them. The function H was then solving a discretized Dirichlet problem. 

In particular, the use of H made the estimation of F on the free arc possible. More 
precisely, F was related to the square root of modified harmonic measures (see Propo- 



sition 3.2 of [DCHNTni). This fact was crucial in the proof of Theorem [2I] |D(]HNin) . 
Without entering into details, let us say that in Dobrushin "domains" {G,a,b), the prob- 
ability at criticality for a site x on the free arc to be connected to the wired arc (which is 



the modulus of the observable, thanks to Lemma 2.2) is on the order of the square-root 
of the harmonic measure of the wired arc seen from x. Equivalently, for a dual site u on 
the wired arc, the probability of being dual-connected to the free arc is of order of the 
square-root of the harmonic measure of the free arc seen from u. We refer to |DCHN10| 
for additional details on these facts. 

We will be using this fact for two nice infinite Dobrushin domains: 

• The infinite strip = Z x [0, n] of height n. Denote by (pg°°p°^ the random-cluster 
measure with parameter p, free boundary conditions on the bottom and wired bound- 
ary conditions on the top. The probability at criticality for a dual-site on the top to 
be dual-connected to the free arc is of order l/\/n (since the harmonic measure of 
the free arc is 1/n, via the Gambler's ruin). 

• The upper half-plane H. Denote by (j>^°^ the random-cluster measure with parameter 
p, free boundary conditions on Z_|_ = {0, 1, • • • } and wired boundary conditions on 
Z_ = {• • • , —2, —1,0}. The probability at criticality for the dual site adjacent to — n 
to be dual-connected to the free arc is of order 1 / ^/n for the same reason as for the 
strip. 

When p < Pc, the fermionic observable Fp can be defined in these two domains even 
though they are infinite (see [BDClOb] ). In the strip, 7 goes from —00 to 00, while in 
H, it goes from to 00. One should be careful about the definition of the winding since 
Cf, does not make sense: the winding is fixed to be equal to on edges of the free arc 
pointing north-east. Since the observable in infinite volume is the limit of observables in 
finite volume, it still satisfies the properties of the previous section. 
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Figure 2.3: Left: The strip. Right: The upper half-plane. In this case, r is the hitting 
time of grey edges. 



2.3 Integrability relations of the fermionic observable away from criti- 
cality 

Away from the critical point, the primitive H is not available anymore. Nevertheless, Fp is 
massive harmonic inside the domain. In fact, Fp satisfies very explicit relations on the free 
arc of the domain, as well. (On the wired arc, such relations are not available.) Indeed, 
Case 3 of Lemma 4.4 of |BDC10bj says that 

cos 2a r . N / M cos(7r/4 + a) , , , . 

^^^^(--)= 2(l + cos(./4-a)) [^^(-")+^^(-"^] + l + coi/A-a) ^^^'-^ " ^^^^^^ 

= (2.9) 

for X on the free arc (excluding in the case of the upper half-plane). When p = pc, 
the sum of the coefficients on the right equals 0, which means that the observable has an 
interpretation in terms of reflected random walks. This relates to the discretization of the 
Riemann-Hilbert boundary problem, and it provides an alternative strategy to handle the 
scaling limit of the observable. 

Away from criticality, we can also interpret these relations in terms of a random process. 
Define the Markov process with generator Ap, which one can interpret as the random walk 

of a massive particle. We write this process {Xn\mn^) where is a random walk with 
jump probabilities defined in terms of Ap — the proportionality between jump probabilities 

(p) 

is the same as the proportionality between coefficients — and m„ is the mass associated 
to this random walk. The law of the random walk starting at an edge x is denoted Pp. 

In order to simplify the notation, we drop the dependency in p in {Xn \'m^n^) and simply 
write {Xn,mn). Note that the mass of the walk decays by a factor cos 2a at each step 
inside the domain, and by some other factor on the free arc. 

Define r to be the hitting time of the wired arc, more precisely, of the set d of medial 
edges pointing north-east and having one end-point on the wired arc (the grey edges in 



Fig. 2.3). The fact that ApFp = for every edge x ^ d implies for any t > that 

Fpix) = E^[Fp(XMr)mMr]. (2.10) 



Since rrir < 1, Fp{Xt/\T)'mt/\T is uniformly integrable and (2.10) can be improved to 

Fp{x)=E;[Fp{Xr)mr]. (2.11) 
This will be the principal tool in our study. 



Proposition 2.5. Let A > 0. There exists Ci = Ci(A) such that for every n > and 

'A5rr(o^^^+^) ^ %■ (2.12) 



Pc>P>Vc-^ 
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There exists C2 = C2 > such that for every n > and p < pc — 

^Zp^iO^zn + Z) < (2.13) 

Proof. In both cases, we study the probabihty for a point on the free arc to be connected 
to the wired arc. In particular, Lemma 2.2 imphes that the quantities on the left of (2.12) 
and (2.13) are equal to |-F(eo)| (or -F(eo) in this case, since the winding is fixed on the 



boundary). Moreover, (2.11) allows us to write 



Let us first deal with (2.12). Note that Fp(Xt-) = Fp{in) = Fp*(eQ), where p* is the 
dual parameter. Hence 



°^'°°(Oom + Z) 



E°[Fp*(eo)m, 



■"pL-' P 

-00,00 
,P 



using Lemma 2.2 again. Since p > pc 



cos 2a 



m 



ip) 



is larger than 1 — c{X/n)^ for 



some c > 0. Using the upper tail of r/n?, we deduce that 

E''[mr] > C 
for some C = C(A). In addition to this, 

C 



Sn,P* 



n 



where we used the standard estimate of the probability at criticality, together with p* > pc- 
The two inequalities together yield (2.12). 

Let us now turn to (2.13). When p < pc — C\/logn/n, we use the expansion of a near 
Pc and cos 2a < 1 — c(logn)/n^ (for some constant c = c(C)) to deduce 



n 



for c' = c'(C). In order to conclude, c' can be chosen larger than 4 by tuning C. 



□ 



The previous proof of (2.12 ) was based on a comparison with the estimates at criticality: 
when p > Pc — X/n the connection probabilities are of the same order as the critical ones. 
We push this reasoning further in the following proposition. 

Proposition 2.6. For any A > 0, there exists C3 = C3(A) > such that for every n > 
and p> Pc - ^, 

</.^^(noZ_) > ^ 



(2.14) 



Let us first prove an easy yet slightly technical result. It should be compared to 
Lemma 12.21 



Lemma 2.7. Let u be a dual vertex adjacent to the wired arc ofM, 

Fp{eu) ^ 0H^('u A Z+), 

where Cu is the edge pointing north-east and adjacent to u, and x means that the ratio is 
uniformly bounded away from and 00. 
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Proof. If V is the vertex of the medial lattice on the left of tt, relation (2.5) around v 
gives F{NW) + F{SE) = e'-°'{F{NE) + F{SW)), where edges are indexed with respect to 
the direction they are pointing to (see Fig. 2.3). Since we know the complex argument 
modulo vr of the observable, we can project the relation on e~*'^''^M, to find 

e'^/^F{NW) - cos(^/4 - a)iF{SW) = cos(7r/4 + a) F{NE). 

Now, the argument of the observable at NW and SW is in fact determined, since the 
winding on the boundary is deterministic (it equals — 7r/2 for NW ^ and — vr for SW). 
Therefore, Lemma [2 . 2| implies 

eW4^(W) = \F{NW)\ = <;^(nAZ+) 

iF{SW) = \F{SW)\ = 0^;^(n-l 

So, using the fact that NE = e^, we get 

cos(7r/4 + a)F(e„) = 05|j^(n A Z+) - cos(7r/4 - a)(/)^^(M - 1 A Z4.). 

Now, (t)^^{u — 1 A Z+) < <l)^^{u A Z+) thanks to the comparison between boundary 
conditions. We deduce 

1 - cos(7r/4 - a) 0,00. A7/\^t?(\^ ^ aO,oo/ ^ 

cos(7r/4 + Q) '^'P ^ ^' ~ ^ - cos(7r/4 + Q) "^'J'^ ^' 

which is the claim. □ 

We are now in a position to prove the proposition. 



Proof of Proposition 2.6 Fix n > and p > Pc — - and denote the fermionic observable 



in (H, 0, 00) by Fp. Lemma 2.7 implies 



Fp{n) = E;[Fp{Xr)mr] x E^i^P^^^iXr A Z+)m,]. (2.15) 



We know that 



6^^;(nAZ+)>C3/ 



hence (2.15) implies 

Fp{n) > C^E;[(t>''^'^{XrAZ+)mr] > C^Cs E;[ mPJ y^l] , (2.16) 

with two universal constants C3, C4 > 0. 

Therefore, it is sufficient to prove that is not larger than n and that rrir is larger 
than some constant e with probability bounded away from uniformly in n. The second 
condition can be replaced by the event r < for instance. 

The walk Xt away from the real axis is just simple random walk, while on the free arc 
it has some outwards drift. So, it is sufficient to prove that Xt exits [0, 2n] x [0, n] through 
[0, 2n] X {n} in fewer than n^/2 steps with probability larger than some constant c > not 
depending on n. Indeed, it has then a uniformly positive probability to exit the domain in 
fewer than n^/2 additional steps and to satisfy \Xt-\ < n. 

Consider (^t)f<n2/2 = (^ti ^i)t<n2/2 conditioned on the event that (Xt)j<„2/2 visits 
the free arc fewer than n times. The probability that the first coordinate is less than n for 
every t < n^/2 is bounded away from uniformly in n (since the number of visits of (Xt) 
to the free arc is less than n, (At) can be compared to a symmetric random walk with a 
deterministic drift of order rn for r < 1). Now, conditioned on the visits of (Xt) to the free 
arc, {At) and (Bt) are independent. Thus, (Bt) is a random walk reflected at the origin 
conditioned on the fact that it does not visit more than n times. In time n^/2, it reaches 
height n with probability bounded away from 0, uniformly in n. The claim follows. □ 
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2.4 Proof of Theorem [Ql 



We first prove crossing probabilities in rectangles with specific boundary conditions. Then, 
we use these crossings to construct crossings in arbitrary rectangles with free boundary 
conditions. 

Crossing in rectangles with Dobrushin boundary conditions. Let us first use the 
estimates obtained in the previous subsection to prove crossing probabilities in the strip 
and the half-plane. The proof follows a second moment argument. 

Proposition 2.8. Fix A > 0. There exists Cg = Cq{X) > such that for every n > and 
P>Pc- ^, 



and 



(l>Zp^{[-n,n]^in + Z)>C^ 



'/'Hp'"([3^'4n]oZ_)>C6. 





A 



3n 4n 



Figure 2.4: The two crossing events of Proposition 2.8 



Proof. We present the proof for Sn (a similar argument works for H). Let be the 
(random) number of sites on [—n, n] which are connected by an open path to in + Z. 
Proposition |2.5| implies that 

'A5rr(^) = E '^5rr(^^^^+^) ^ (^n+i)^ > 2c,v^. (2.1?) 

x£[—n,n] 

Moreover, for p < Pc, 

The right hand side is a quantity at the critical point and was already studied in the proof 
of the main theorem of |DCHN10] (in fact, only closely related quantities were studied, 
but the generalization is straightforward). In particular, it was proved in that article that 

Cauchy-Schwarz thus implies that 

'A5rr([-^' n] o m + Z) = c^s^^f^iN > 0) > ACf/C, , 
uniformly in n. For p > Pc the result follows from monotonicity. □ 
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It is now easy to reduce crossing probabilities in the strip and the half-plane to crossing 
probabilities in (possibly very large) rectangles. The idea is that a crossing cannot explore 
too much of the strip or the half-plane, since there exist slightly supercritical dual crossings 
preventing it. 

Proposition 2.9. Fix A > 0. There exist C7 > and M > such that for every n > 
and p> Pc- ^, 

€~MS:l"Sk.(^--^ n] o .n + Z) > 

and 

Proof. As before, we do this in the case of the strip. Fix M large enough so that, 
at criticality, the probability that there exists a vertical dual crossing with free boundary 



conditions of [n, Mn] x [0, n] exceeds 1 — Cg/S (use Theorem 2.1 to prove this fact). Then, 
with probability Cq/3, there will exist a crossing of [—n,n] to in + Z and two dual vertical 
crossings in [n, Mn] x [0, n] and [—Mn, —n] x [0, n]. The domain Markov property and the 
comparison between boundary conditions imply the result. □ 



Crossing in rectangles with free boundary conditions. A consequence of Proposi- 
tion 2.9 is the existence of crossings inside a box with free boundary conditions everywhere. 



Indeed, although the previous result only deals a priori with domains where a part of the 
boundary is already wired, this condition can be removed. 



Proposition 2.10. Fix A > 0. There exist Cs,M > such that for every n > and 
'^['-Afn,Mn]x[o,n],p([~^'^'^"'] ^ [O'"'/^] crosscd Vertically) > Cg. 



P>Pc-^ 



Proof. Fix M so that Proposition 2.9 holds true. Let An be the event that [—Mn, Mn] x 
[0,n/2] is crossed vertically. We have for every n > 0, 

,{i-M)n,{i+M)n (a \ ^ r 
^[-Mn,Mn]x[0,n],py^'"-) - ' 



Let Bn be the event that [—Mn,Mn\ x [n/2,n] is dual-crossed horizontally. Theorem 2.1 
implies that 

.{i-M)n,{i+M)n , . . 

^[-Mn,Mn\x[0,n\,p\^'^\^'^) - ^ 

for some constant c > 0, uniformly in n and p < pc- Now, 

,{i~M)n,{i+M)n 



V'[-Mn,Mn\x[0,n],p\-^n) ^ "/^[-Mn.Mn] X [0,n],p*^^" I 

,(i~M)n,{i+M)n „ ^ 

- 'i^[-Mn,Mn]x[0,n],pV^" ' ' 

_ Ai-M)n,{i+M)n /r i . n , (i-M)n,{i+Af)n ^ 
~ '^[-Mn,Mn]x[0,n],pV-Dn|^n; • 'P[-Mn,Mn\x[0,n],py^^l 

> c-Cj. □ 
We now prove that crossings of rectangles of any aspect ratio also exist. 



Lemma 2.11. Fix A > and k > 0. There exists Cg = Cg{K, A) > such that for every 
X 
n' 

f-n,(K+i)n]x[o,n],p([0>'^^] >< [0' '^l crosscd horizontally) > Cg. 



n andp > Pc- ^, 
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n 



Men 



event A 



event Bh 




{k + l)n 



ekn e(k + l)n 



Figure 2.5: The intersections of the events A and create a crossing of the rectangle 
[— n, n] X [0, Kn]. 



Proof. Fix M = M(A) as in Propositions [Z9| and [2lo| Let e = 1/(2M)2. Let A be 
the event that there exists a crossing from [—en, en] to zMen + Z, and let Bk be the event 
that there exists a path in Z x [0, Men] from [(/c + l)en, [k + 2)en] to [{k — l)en, /cen]. See 



Figure 2.5 We have 



0f_„,(«+i)„]x[o,n],p([-"''^] X [O''^"] is crossed horizontally) 

k/e— 1 



> 



-n,(K+l)n]x[0,n],p I 1^ fl 
\ fc=0 

ft/e— 1 

,(^) J] </'f_„,(.+iHx[o,n],p(^fcl^'^r,r</c) 



-n,(K+l)n] X [0,n] ,pV 



fc=0 



Furthermore, 



^[— ri,(At+l)n] X [0,n] ,p -n,n]x[0,n/(2M)],p 

Now, since Me = 1/(4M), the event A in [— n, n] x [0, n/ (2M)] corresponds to the existence 
of a crossing from the bottom to the middle, but with the additional constraint that it 
starts between [—en, en]. A union bound, comparison between boundary conditions, and 



Proposition 2.10 imply that 

'^[-n,n]x[0,n/(2M)],p(^) > e'^[-„,n]x [0,n/(2M)],p([-^' X [0, »^/(4M)]) > e • 

Furthermore, 

4>[-n,n]x[-n,{K+l)n],piBk\A,Br,r < k) > (l)'([k-M)en,(k+M)en]x[0,M£n],p(^k) > Cj , 

using the comparison between boundary conditions and the half-plane case of Proposi- 



tion 2.9 Altogether, we obtain that 



^'P_,,,(«,+i)„]x[o,n],p([-"''^] X [0''^^] is crossed horizontally) > CsCy^^, 



and the lemma is proved. 



□ 
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event Ctop 




event Cbottom 



Figure 2.6: The five events involved in the proof of Theorem 1.3 



Proof of Theorem 

[en, (1 — e)n] x [en, 2en] and [en, (1 



and Atop be the events that 



Fix e < 1/(4M). Let Abottom 

e)n] X [{p — 2e)n, {p — e)n\ are crossed horizontally. 
Let B be the event that [en, (1 — e)n] x [en, (p — e)n] is crossed vertically. Let Chottom 
and Ctop be the events that [en, (1 — e)n] x [0, 2en] and [en, (1 — e)n] x [(/> — 2e)n, pn] are 
crossed vertically. See Figure 2.6 By Lemma 2.11[ the events Abottom, Atop and B have 
probability bounded away from uniformly in n. The FKG inequality implies that their 
intersection also has this property. Now, conditionally on Abottom, the event Cbottom has 
probability larger than the probability that there exists a crossing in [en, {l — e)n] x [0, 2en] 
with wired boundary condition on the top and free boundary condition on the bottom. 
Proposition 2.9 implies that this probability is larger than C7 since (1 — 2e)/(2e) > 2M 
(the important thing is that the rectangle [en, (1 — e)n] x [0, 2en] is wide enough). The 
same reasoning can be applied to Ctop, ergo the claim follows. □ 



2.5 Proofs of Theorems O and fOl 

First of all, a standard reasoning described in Step 2 of the proof of Theorem 1 of [BDClOb] 
shows that equation (2.13) implies the following lemma: 

Lemma 2.12. There exists Cio > such that 

(pp{0 ^ d[-n,n]'^) < Cwn'^ (2.18) 
for every n large enough and every p < Pc — C'lo -. □ 
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Proof of Theorem \1.2\ Fix Cio > as defined in Lemma 2.12 Theorem 1.3 imphes 



the lower bound trivially. For the upper bound, it suffices to show that for any k > 0, 

<A[_„,„]x[-«nAn],p([-'^''^] ^ [-i^n,nn] is crossed horizontally) ^ 

whenever {n,p) — )• (oo, 0) with p < Pc — Cio\/logn/n. Fix e > and k > 0. 

Take some 5 > to be fixed later. Let be the event that [—(1 — 6)n, (1 — 5)n] x 
[—KTi, (k — 26)n\ is crossed horizontally, and A^°^^°^ be the event that [—(1 — 5)n, (1 — 
5)n] X [—{k — 2(5)n, Kn] is crossed horizontally. Furthermore, let be the event that 
[— (1 — 5)n, (1 — 5)n] X [— (k — (^)n, (k — contains a cluster of diameter 5n. Notice that if 
the rectangle [—77 - , n] x [— /tn, Kn] is crossed horizontally, then A'n^ or or Bn occurs. 



Theorem 



2.1 



implies the existence of > such that the probability of A^n^ (and 
similarly for yl^*"'") with wired boundary conditions is smaller than e/3 for any p < psd 
and n > 0. This will be our 5. 

Define Cn to be the event that the annulus 

Sn '■= [— n, n] X [— Kn, nn] \ [—(1 — 5)n, (1 — Sn)] x [— (k — 5)n, {k — 6)n] 

contains a closed circuit surrounding the inner box. Note that there exists r] > such that 



thanks to Theorem |2 . 1 1 again . Since C„ is decreasing and Bn depends only on edges inside 
[— (1 — 5)n, (1 — 6)n] x [— (k — S)n, {k — 5)n], we obtain 



Therefore, 



^l,[-n,n]x[-Kn,Kn]i'^n\Bn) > </'p_5^(Cn) > V- 



^p,[-n,n]x[-Kn,Kn]i-^n) < 4'p^[-n,n]x[~Kn,Kn]i-^ri r) Gn) 
— ^v.\-n.n]x\-Kn.Kn](^n\Cn) 



p,[—n,n] X [—Kn,Kn] " 
— ^p,[-n,n]x[-Kn,Kn]i'^n)- 



By a union bound. Lemma 2.12 and the definition of Cio imply that 



<[-n,n]x[-Kn,.n](^") when n ^ 0. (2.19) 

Therefore, </>J^[_^„],[_,„^,„](i?„)^0. 

Summarizing, each of A^^^, A^^^°"^ and Bn has probability less than e/3 for n large 

enough, which concludes the proof. □ 



Let us now turn to the proof of Theorem 1.4 We have just proved that, for p > and 



e > 0, there exists c = c{e, p) such that for any n > jj-^wlog ^ 



'^p,[0,n]x[0,pn](C?»([0,n] X [0,H)j < £• 

The next lemma asserts that crossing probabilities in fact converge to very quickly as 
soon as n is larger than the correlation length. 
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Lemma 2.13. For any p < pc, there exists L{p) such that 

1 



Pc-P 



< L{p) < 



log 



1 



c{Pc -P)\ Pc-P 



and 



_2k 



<[0,2'=L(p)]x[0,2^+iL(p)](C/.([0,2^^(p)] X [0,2^'+iL(p)])) < 
for any k > 0. 
Proof. For n > 0, let 

Un :=max|(/.^ [o „]x[o 2n](Ch([0,n] x [0,2n])) , <^J,[o,n]2 (C/i([0, n]^)) | . 

We are going to show that 

U2n<'25ul. (2.20) 

Fh'st, cutting vertically the domain [0,2n]^ into two rectangles, together with compar- 
ison between boundary conditions, implies that 

<[o,2n]2(C/.([0,2n]2)) < ^1 [Q„j^[Q2n](Ch([0,n] X [0,2n]))' < ul . (2.21) 

Second, cutting vertically the domain [0, 2n] x [0, 4n] into two, together with comparison 
between boundary conditions again, implies that 

<[0,2n]x[0,4n](C/^([0,2n] X [0,4n])) < <Ap,[0,n]x[0,4n](Ch([0,n] X [0,4n])'' 
Now, consider the rectangles 





= [0,n] 


X 


[0, 2n] 


R2 


= [0,n] 


X 


[n, 3n] 




= [0,n] 


X 


[2n, 4n] 


R4 


= [0,n] 


X 


[n, 2n] 


R5 


= [0,n] 


X 


[2n, 3n] 



These rectangles have the property that whenever [0, n] x [0, 4n] is crossed horizontally, at 
least one of the rectangles Ri is crossed (in the horizontal direction for R2 and R3, 
and vertically otherwise). We deduce, using the comparison between boundary conditions, 
that 

<^J,[o,n]x[o,4n](Ch([0,n] X [0,4n])) < 5Un, 

and hence 

<[o,2n]x[o,4n]fe([0,2n] X [0,4n])) < {5unf. (2.22) 



Combining (2.21) and (2.22), we obtain (2.20). Iterating that, we easily obtain that, 
for every A; > 0, 



25 U2kn < (25 Ur, 



1.2 



By Theorem 



then Un satishes 



ifp < pcandn > ^ Jlog^, where c = min{c(l/100, 2), c(l/100, 1)}, 



25 n„ < 1/e . 



Therefore, the lemma follows for L(p) = — — a /log ^ . 

' Pc-py pc-p 



□ 
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Proof of Theorem 1.4[ Fix p > pc- Let 

Rk := [Q,L{p)2^] X [-L(p), L(p)(2^+^ - 1)] if A; is even, 

and 

Rk := [0,L(p)2'=+i] X [-L{p),L{p){2^ - I)] if it is odd. 

Define Ek to be the event that Rk is crossed in the long direction. The FKG inequality 
implies that 



bliS) > </>°(0 o {L{p)} X [-L{p),L{p)\) ■ W ct>l{Ek) 

k>0 

> ^^cl^l{0^d[-L{p),L{p)]').ll{l-e 

> c{L{p))-''\ 



k>0 



where c > 0. We used Lemma |2.13| to get the second line, and the lower bound of Theo- 

□ 



rem 



1.2 and (1.1) to get the third inequality. 



3 Near-critical behavior: a fascinating self-organized near- 
criticality emerges 

As promised, we start the discussion with the near-critical regime in standard percolation. 
Our goal here is to provide a self-contained explanation of the fact that the near-critical 
correlation length and the behavior of the percolation probability 0{p) as p \^ pc are 
governed by the number of pivotals at criticality. Then we explain why this picture must 
be flawed in the case of random-cluster models. 



3.1 Pivotal points govern the near-critical regime of percolation 

Recall that the correlation length is defined as follows: fix some small e > 0. For any n > 1, 
let Rn be the [0, n] x [0, n] square for bond percolation on Z^, and for any p = pc + Ap > pc, 
define 

Le{p) = L{p) := inf {n > 1 s.t. Pp(C^(i?„)) > 1 - e} , 

where Pp is the probability measure of bond-percolation with parameter p. Let us start by 
explaining the fact that things look supercritical above L{p). For n > L{p), the probability 
to have a left-to-right crossing in the box [0,n]'^ is larger than 1 — e by definition. Russo- 
Seymour- Welsh theory (see |Gri99| Lemma 11.73]) implies that the probability to have a 
left-to-right crossing in the rectangle [0, 3n] x [0, n] is greater than 1 — g{e), where g{e) 
goes to as e — )■ 0. Then, using well-known arguments, one can show that the "geometry" 
of ojp above L(p) stochastically dominates a certain supercritical percolation model of 
parameter 1 — (j)(g{e)), where (j){x) goes to as x — )• 0. In this sense, things indeed look 
supercritical above L{p). This step would have been simpler if one had worked directly 
with long rectangles [0,3n] x [0, n] in the definition of L{p) instead of the symmetric Rn- 
However, the symmetry of i?„ will be relevant to the explanation of the following opposite 
fact: things look critical below L{p). For this, we need to obtain RSW estimates for the 
dual percolation. By planar duality, if n < L{p), then the probability to have a top-to- 
bottom dual crossing in R^ is greater than e. Russo-Seymour- Welsh theory then implies 
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that the probabihty to have a dual left-to-right crossing in [0, 3n] x [0, n] is greater than 
V'(e) (where ip{e) — t- when e — ?■ 0). This means that at scales smaller than L{p) we do 
have the uniform crossing probabilities that make the configuration look critical. 

Goal: How can one estimate L{p) as a function of \p — Pc\'^ 

As explained in the Introduction, it is quite simple to guess what L(p) should be. 
Indeed, the quantity is naturally related to thermodynamical quantities of the model via 
Kesten's scaling relations: 

Theorem 3.1 ( |Kes87| . see also |Wer07| fNolOS] for modern expositions, and |GPS| for an 
alternative approach). For L{p) = L^lp), one has 

L{pfai{L{p)) X (3.1) 

\P-Pc\ 

e{p) X ai(L(p)), (3.2) 
where the constants in x depend on e > 0. 

Sketch of proof. Let us start with the first relation. In the following, let denote the 
indicator function of the left-to-right crossing event in the domain Rn- 

The intuition was already given in the Introduction: for crossing events of scale n (i.e., 
for rectangles of diameter x n), there are Q{n?a4{n)) points on average which are pivotal 
for the crossing event. Now, in the standard monotone coupling, from ujp^ to LVp, each of 
these pivotal points flips with probability of order |p — Pc|- Therefore, it is tempting to 
believe that as far as n?a4{n)\p — Pc\ ^ 1, it is unlikely for the crossing event to change, 
while once n'^a4{n)\p — Pc\ ^» 1, many pivotal points are flipped and things should start 
being highly connected. 

A few things are a bit "fishy" in this intuition: one of them is that one might have 
/n(wp) = 1 and fni^^pc) = together with the fact that from ujp^ to ujp, none of the initial 
pivotal points for /„ had been switched: the pivotal switch might happen at a point that 
was not pivotal originally. On the way from Up^ to ujp, if one stayed at equilibrium, as 
it is the case, for example, in dynamical percolation, one would still be able to conclude 
something based on such considerations, but one difficulty here is that as we follow the 
monotone coupling, we leave critical percolation. 

The nice idea from |Kes87] to overcome this near-critical bias is to apply Russo's 
formula simultaneously to the crossing event as well as to the four-arm event. Indeed, for 
the crossing event, one can check that as long as n < L(j)) = L^{p), 

^^p[fn] = Yl IPp[2; is pivotal] 

sites X 

X n'^a^{n) . (3.3) 

To go from the first line to the second one is in fact non-trivial: one needs to prove that 
below the correlation length, the main contribution in Russo's formula comes from bulk 
points rather than boundary points. The technology involved here is quasi-multiplicativity 
and separation of arms. One can prove these even if we are not at the critical point, since 
as far as n < L(p), one still has RSW estimates both for the primal and the dual model. 
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Now the key observation is the following one: 
d 



dp 



al{n) 



< 



^^Pp[2; is pivotal for the 4-arm event] 



using Russo's 
formula again 



< 0(1) a\{n/2)) ^p[x has the 4-arm event to distance n/3] 

|a:|>2ra/3 

<0{l)al{n)n^al{n). (3.4) 

The second inequality uses quasi-multiplicativity along with a dyadic summation to show 
that the main contribution arises from large-scale pivotal points. The third inequality uses 
quasi-multiplicativity again. 

The combination of (3.3) and (3.4) implies that for all n < L{p), the variation of 
p I—;- log(Q4(n)) is controlled (up to constants which depend only on e) by the variation of 
p I— )• Ep [/„] which is of course bounded. This implies that 

a^{n) X 04(12) , (3.5) 

where the constants involved in x depend only on e > 0. 

Now integrating (3.3) and using (3.5), one can conclude about the correlation length: 
for n < L{p), 



X |p — pc\n'^aii{n) . 



In particular, for n = L{p), since [/„] — Ep^ [/„] x 1, we obtain our desired estimate (3.1 ). 

□ 

The philosophy behind the previous argument is that, as far as few pivotal points are 
touched, percolation remains "critical". This is what we called near-critical stability in 
Phenomenon 1.5 In particular, critical exponents remain unchanged below L(p). We have 
seen this in (3.5) for the case of 04(11), but the same argument works for the one-arm event 
Qi(n): namely, for n < L{p), 

(n) X ai(n) . 

This immediately implies the second scaling relation: 
e{p) :=Pp[Oooo] 

^Fp[0^dB{0,L{p))] 



since above L{p), there 
is a "dense" infinite cluster 



amp)) 

ai{L{p)) . 



□ 



The knowledge of critical arm-exponents ^i(l) and (,4(1) for site-percolation on the 
triangular lattice T (see |SWOH rLSW02j ) allows for an estimation of L{p) and 6{p): 

Lj{p) = (p-p,)-4/3+o(l) 
0j{p) = (p-p,)5/36+o(l) 



as p\pc = 1/2. 
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As we discussed at length in the Introduction, the first scahng relation (3.1) does 



not hold in the FK-Ising case (see (1.4) and Phenomenon 1.6). Nevertheless, note that 
the heuristics of the second scaling relation should work in great generality. Contrary to 
the first scaling relation, it gives the right prediction for the correlation length. Indeed, 
the thermodynamical quantities L{p) and 9{p) have their analogues in the FK-Ising case. 
Onsager's determination of the magnetization, together with the Edwards-Sokal coupling 
implies that 

el¥{p) := <Pp,2{0 ^ oo) X \p-pc\^/^ (3.6) 

and 

^p^M^ ^ d[-n,nf) X (3.7) 

From these two relations, the second scaling relation (which does not harness any pivotal 
event) implies that the correlation length should behave as l/\p — Pc\ for FK-Ising, which 



is the right prediction. Also note that (3.7) has been proved using conformal invariance 



techniques in |CHI12j . It would be interesting to make sense of the second scaling relation 
in the FK-Ising case in order to provide a derivation of the exponent 1/8 for the magneti- 
zation which would be independent of Onsager's computation. Half of this is achieved by 
Theorem 11.41 

Let us conclude this paragraph by mentioning another result proved in |Kes87| . The 
correlation length is sometimes defined as the inverse rate of exponential decay of the 
connectivity probabilities in the subcritical regime: 

lim log Pp [(0,0) o (n,0)l. 



L{p) n-5-oo n 

Kesten proved that this definition is consistent with the definition in terms of crossing 
probabilities. Nonetheless, this correspondence is known only in the percolation case and 
is not fully understood in the FK-Ising case. Furthermore, this definition is less tractable 
when studying near-critical regimes. It is therefore useless to consider existing estimates 
on the inverse rate of exponential decay for the FK-Ising model (see |BDC10b) e.g.) to 
prove that the correlation length defined as before behaves as \p — Pc\~^ ■ 

3.2 Specific heat of the random-cluster model 

We would now like to understand the behavior of random-cluster models when p varies 
from to 1. The first non-trivial effect which occurs in the near-critical random-cluster 
model is the fact that the derivative of the edge-intensity blows up around pc for q > 2. 
This implies that edges appear much faster in any monotone coupling near pc than they 
do for percolation. Define the edge-intensity function as follows: for all p G [0, 1], let 

^{p) ■= (l>p,q{e is open) , 



where e is any edge of Z^. It is not hard to check that at the critical (and self-dual) point 
Pciq) = one has 

A quantity relevant to us is the derivative in p of the edge-density dX{p)/dp. It corre- 
sponds to the average rate at which new edges appear in any possible monotone coupling 
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i^^p,q)p£[o^i]- For integer q > 2, this quantity turns out to be linked to the so-called specific 
heat of the (7-Potts model, defined by 

^Potts(q)(^) := -f]^^E^[H] = /32Var/3[/7] , (3.8) 

where H is the total energy of the system, and where the second identity can be easily 
proved by manipulating the partition function (recall here, e.g., that IE/3[-ff] = 

— (9^ log By the Edwards-Sokal coupling, the partition function equals 

the partition function Zp^q of the FK{p, q) model with e~^^ = 1—p, hence the specific heat 
can also be defined for the FK{p,q) model, to be denoted by C^^''{p), using the second 
derivative of the free energy log Zp^g. This works now for any q > 0, except that the q = 1 
percolation case (where Zp^i = 1) needs a limiting procedure q ^ 1, which yields the 
second derivative of the expected number of clusters per site, or £Ep[l/\Ca;\]/dp^, where 
Cx is the cluster of any fixed site x. 

The relationship between the derivative of the edge-intensity and the specific heat per 
site is detailed in |GH| . but the point is that for q > 2, they are within bounded factors 
from each other. So, for q = 2, from the results on the specific heat of the Ising model 
known since |Ons441 IFF69] . one obtains in the infinite volume case that, at p = Pc(2), 



d 
dp 



P=Pc 



X{p) = 00 



More precisely, one can extract the following logarithmic behavior around pc- 

-T{p) ~ a log ■ 



dp IP - Pel 

as p — Pc- The finite- volume study of the specific heat ( |Ons441 IFF69| ) leads to the 
following estimate: let T„ be the torus 1? jnl? and let p 1— )• Inip) denote the edge-intensity 
for random-cluster model on T„, then 



dp 



In{p) X logn. 



P=Pc 



The extension to planar domains il„ := ^Z^ n Q, will be carried out in |,GHj . based on the 
recent results from |HonlOj and [ BgTOSl [BHT09| . 

In conclusion, these estimates show that in a window of size n, as one raises p near pc, 
more edges will suddenly arrive. Nevertheless, the discrepancy is only logarithmic, hence 



it is not sufficient to explain the error (1.4) in the power law of the correlation length. 



3.3 The monotone increasing Markov process on cluster configurations 

Let us present a monotone coupling /i between random-cluster models with fixed q > 1, 
which was first considered by Grimmett in jGri95| . It turns out that it is non-trivial to con- 
struct explicitly such a measure ^ (note that in contrast the existence of abstract monotone 
couplings follows easily from a generalized Strassen's theorem). Instead of constructing ex- 
plicitly a coupling fi, Grimmett obtains it for a graph G = {V, E) as the invariant measure 
^ of a natural Markov process on the space VL := [0, l]'^. This technique is usually employed 
to simulate a single Gibbs measure and we propose to provide this example first. 
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Heat-bath dynamics. Assume we wish to simulate the random-cluster measure (j^cpq 
on the graph G = {V, E). First note that the random-cluster model with parameters [p, q) 
has the following property: 



Property 3.2. For any edge e = [xy], 
^%,p,q{^ is open \ uj on G \ {e}) 



p , if X y in G\{e} 
, ,? — X- , otherwise . 



This "almost local" rule for the conditional law of an edge e knowing the external 
environment enables us to consider the heat-bath dynamics: 

Definition 3.3 (Heat-bath dynamics or Sweeny algorithm). Let G = {V,E) he any 
finite graph. The random- cluster heat-hath dynamics on G is defined as follows: each edge 
e € E is updated at rate 1 (the exponential clocks heing independent) and when a clock rings 



at e, its status ijj{e) is resampled according to the conditional law given in Property 3.2 



It is straightforward to check that this dynamics has the following properties: it is a 
reversible Markov chain with state space {0, 1}^ and its invariant measure is 4>Gpq- This 
dynamics has been studied both for theoretical reasons (see |Gri95| IGri06| ) and practical 
ones (see |DGS07| for a good account of recent works). For instance, via the Edwards-Sokal 
coupling, it turns out to provide a faster way than classical Glauber dynamics to sample 
Ising models (at least in dimension d = 2). More sophisticated algorithms are known 
for integer values of q (for example, the Swendsen-Wang algorithm). However, the above 
dynamics has the advantage that it works for all real values of q and is probably more 
tractable for rigorous analysis. 

Grimmett's dynamics and monotone coupUng. We now describe briefly Grimmett's 
monotone coupling (see |Gri95[ [Gri06] for a detailed exposition). Let G = (V, E) be a finite 
subgraph of I? and Q be the space [0, 1]^. Each Z G Q decomposes into a monotone family 
of edge configurations (ti;p)o<p<i = (a;p(Z'))o<p<i, where for each p £ [0, 1] and any e G E: 

ujp{Z){e) := lz{e)<p- 

The goal is to find a measure fi = ^iq on in a such a way that all the "projections" 
ijjp{Z) with Z ^ jjL follow the random-cluster probability measure of parameters {p,q) on 
{0, 1}^ with free boundary conditions. It is not hard to see what this Markov process 
{Zt)t>o should be. Assuming that for all p E [0, 1], the projection Up{Zt) is also a Markov 
process and that its invariant measure is (pQpq, then it is natural to expect ujp{Zt) to be the 
heat-bath dynamics that we defined previously. Indeed, if {Zt)t>o is such that each edge is 
updated at rate one (the exponential clocks on each edge being independent), this means 
that simultaneously for all p S [0, 1], the law of the update at e needs to be compatible 
with Property 3.2 For any e = (x, y) G E, let C {0, 1}^ be the event that there is a 



path of open edges in E\ {e} connecting x and y. For any e £ E and any Z G ^1, define 

Te{Z) := inf{p € [0, 1] s.t. ojp{Z) G V^} . 

Assume one is running the dynamics and that at a time t the edge e rings. Let Zt- be 
the current configuration (before the update) . Let Ue be the random variable corresponding 
to the new label at e knowing Zt-. In particular, we know the value of T = Te{Zt-). Since 
p >T is equivalent to ujp{Zt-) G Pg) in order for Zt to match the heat-bath dynamics on 
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the projection ujp{Zt) for all p simultaneously, conditionally on the value T = Te{Zt-), the 
update variable must satisfy 

r lip if p > T , , 

¥[Ue<p] := p (3.9) 

Fortunately, q > 1 implies that this is a valid distribution function, hence we can simply 
define Z^e to be a sample from this distribution. Note that Ue has an absolutely continuous 
part plus a Dirac point mass (for g > 1) on T, namely [T — y^^j^^^^ Jt^r. 

This discussion motivates the introduction of the Markov chain Zt on the state space 
Xa = [0, where labels on the ed ges are updated at rate one according to the above 
conditional law (given hyUg, e £ E(A)). This is precisely the dynamics that was considered 
by Grimmett in |Gri95| (see also |HJL02] where this dynamics was revisited). 

Constructing an infinite- volume version of the previous dynamics is not straightforward. 
Nevertheless, one has the following asymptotic statement from |Gri951 [Gri06j . 

Proposition 3.4 (Infinite Volume Limit |I Gri95| ). For each n>l, let A„ := [— n, n]*^. Let 

^ be some initial configuration in X := [0, ). Consider the random- cluster heat-bath 
dynamics Z^" on A„ with free boundary conditions and which starts from the initial state 
Zq" = ^\A„- Then, as n ^ oo, the process {Z^") weakly converges to a Markov process 
[Zl^'^'^)t>o which starts from the initial configuration Zq'^^ = ^. 

Furthermore, as t ^ oo, Z'^'^^ weakly converges to an invariant measure fj, on X. 

If, in the limiting procedure, one uses wired boundary conditions instead, one obtains 
at the limit a Markov process {ZJ'^'^'^'^)t>o- The processes Z'l'"'^'^ and Zf°'^ might possibly 
have different transition kernels but they both have the same fi as the unique invariant 
measure. 

The underlying dynamics here is non-Fellerian, and the limiting Markov process in 
the above theorem is derived from the monotonicity properties inherent to the dynamics. In 
particular, the relationship between this Markov process and its formal generator (we will 
not write it down explicitly here) would need to be investigated. This seems to be a non- 
trivial task for the present dynamics. Therefore, we will not assume any explicit transition 
rule for the infinite- volume dynamics Z^^'^'^ (or Z^'^^^'^) and will restrict ourselves to the 
"compact case". 



The Markov property of the monotone couphng. Let us prove that Grimmett's 
coupling leads to a monotone increasing Markov process (as p varies) on the cluster 
configurations, i.e., on the space {0,1}^. We will not rely on this Markovian property 
later on. Nevertheless, it provides a nice picture of the self-organization scheme near Pc{q)- 
Namely, as one raises p near pc, new edges arrive in a complicated fashion yet depending 
only on the current configuration LOp. 

Proposition 3.5. Let G = {V,E) be a finite subgraph of 1? . Let Z be sampled according 
to the law /i. Then the monotone family of projections (a;p(Z))o<p<i, seen as a random 
process in the ''time" variable p, is a non- decreasing inhomogeneous Markov process 
on the space {0, 1}^. 

Proof. We wish to prove that conditioned on the projections (a;u(Z))o<u<p, the conditional 
law of the higher configurations (wu(^))p<u<i depends only on a;p(Z). To achieve this, it is 
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enough to prove Lemma 3.6 below. Before stating the lemma, we introduce some notation. 



For p S [0, 1], decompose the configuration Z into the triple [ujp, Z-^, Z^^^) defined as 

ujp = ujp{Za); ^ =\. ,1 • ; ^ 1 n .1 

I 1 otlierwise I U otlierwise 

Note that 

ujp = u:p{Z^P)=uOp{Z>P), (3.10) 
and that Z can be recovered from the triple {ujp^ Z-^, Z'^P). 

Lemma 3.6. Conditioned on the value of the first component Up, the two other components 
Z-P and Z^^ are conditionally independent. 

Proof of the lemma. Fix p G [0, 1] and omit it from the notation uj = uip to make space 
for a time variable t. 

We basically follow the construction of the measure /i as the limiting measure of the 
Markov process Zf, except that we divide the randomness used along the Markov chain 
into three components, the second and third being independent conditionally on the first 
one. Namely, define a Markov process 



K, zp, z>%>o e {0, X [0, X [0, 1] 



E{A) 



where edges are updated at rate one, in such a way that the relations (3.10) between the 
three coordinates hold for all t > 0. To be consistent at t = 0, the process starts either 
from the empty state (wq, Z^^, Zq^) = (0, 1, 1) or the full state (1, 0, 0), where and 1 
denote the vectors all and all 1 respectively. Then, instead of sampling lAf^ directly, let us 
proceed stepwise: first look whether oot^ satisfies or not. If it does, then let ujt{e) := 1 
with probability p. If iVt- ^ "D^, then let (e) := 1 with probability p/ (p + (1 — p)q). This 
is exactly the heat-bath dynamics for (j)Q ^ ^. Note that this part of the dynamics does not 
use at the two components {Z-^, Z^^\ 

Let us describe how to update the component Z^^ . If, after the update, uJti^e) equals 
0, then we fix Zp(e) := 1. Other wise (if Wj(e) = 1), we use the following variable: 

r|^P(Z^P) := inf{n G [0,p] : oj"(Z^p) G 2?e} ■ 

Note that Tf'^{Z^^) = Te{Z) on the event Te{Z) < p. Otherwise {i.e. ujp{Z) ^ Ve), 
we set Tg-^ = p. In either case, it is important here that no information about the third 
component Z''^ has been used. 



Next, recall the update random variable lAe from the previous subsection (see (3.9)). It 
needed as an input the value of T(.{Zt-)- Let lAg"^ be the same random variable here, with 



input the value of T^^{Zf^). Remembering that we are in the case uit{e) = 1, update the 
value of Zf^ as follows, independently of everything: 

zf^{e)^C[Uf^\UfP<p\, 

where C stands for the law of the variable. We define (Z^^) in the same fashion, using 
U^^. In particular, the evolutions of (Zf^) and (Z^^) are sampled out of the evolution 
of (cjf) plus some randomness in each case that are independent of each other, hence the 
conditional independence of Z-^ and Z^p is satisfied. 
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To conclude the proof, one just has to notice that if one defines 



Zf^ if LOt{e) 
else , 



z>p 



then (Zt)t>o is exactly the Markov chain which was considered by Grimmett in |Gri95) . 
(This is not hard to check; an important feature here is that if Tg > p, then the conditional 



law C 



Ue\Ue<p 



does not depend on the exact value of Te, and a similar thing holds 



for lA^^ when Tg < p.). In particular, from |Gri95) . it converges to the unique invariant 
measure /^Aj which inherits its conditional independence property. This finishes the proof 
of Lemma 3.6 and hence of Proposition 3.5 □ 



Proposition 3.7. This Markovian property extends to the infinite volume limit ^ on X = 

Indeed, the same procedure works, but one has to be a bit careful with the initial state 
of our Markov chain: following |Gri95) . with the slightly asymmetric projection convention 
we have chosen, we need to start from the full state (wq, Z^^, Zq^) = (1, 0, 0). In the case 
q = 2, this is not very important since there is a unique infinite volume limit for all values 
of p. However for larger values of q, these considerations do matter. We will not enter in 
more detail here; see |Gri951 IGri06| for a detailed exposition on the infinite volume limit 
of Z together with its projection ojp{Z). 



3.4 Existence of emerging clouds when q = 2 

We now give a concrete manifestation in Grimmett 's coupling of the self-organized be- 
havior appearing in any monotone coupling of random-cluster models, explained in Phe- 
nomenon 1.6 We restrict ourselves to the finite case, since the transition rule for the 
infinite volume Markov process {Zt)t>o has not been established. Let then A be a finite 
box in Z2 (or a torus Z'^/nZ'^). 

To start with, let us define properly the notion of cloud for a finite graph A = (V, E). 
Given a sample Z = Z\ G [0, 1]^ from Grimmett's monotone coupling //a, for an edge 
e G -E, let cloud (e) be the set of edges which appear simultaneously with e: 

cloud(e) ■.= {f s.t. Z{f) = Z{e)] . 

The following proposition gives the first hint of some "non-linear" behavior: 

Proposition 3.8. Fix q> 1. For any N > 1, let ((jjp(ZA)pg[o,i] be a monotone coupling in 
the box A. The probability that clouds of at least N edges appear simultaneously in u}p(Z\) 
at some p G (0, 1) converges to 1 when the size of the box K 1? . 

This proposition is very easy to prove, yet one already sees here that the monotone 
Markovian coupling {oJK,p)o<p<i has a nature that is very different from the q = 1 case. 

Proof. Let us consider the sets Ei^ E2 and E^ in A (which is assumed to be large enough) 
as defined in Fig. |3.1| 

Now let us sample Zq = Z^q according to the invariant measure and let us run 
the dynamics for a unit time. With positive probability, all edges in E are updated and 
their labels at time 1 satisfy the following: all labels in E2 are smaller than 1/4, the edge 
eo = ((0, 0), (1, 0)) gets a label in (1/4, 1/2), and all other labels in Ei U E3 are larger than 
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E, :=U{((0,0,(1,0>} 



1=0 



(horizontal inner edges) 



n-l 



E2:= U {((0.0,(0,^ + 1)), ((1,0,(1,^ + 1))} 



(vertical inner edges) 



Es := {all edges neighboring Ei U E2} \ {Ei U E2) 
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Figure 3.1: The definition of the sets Ei, E2, E^ and the edge eo G Ei. 



3/4. Under such circumstances, all edges e £ Ei \ {eo} are such that Te{Zt=i) = Zi{eo). 
It could be that this situation evolves later on, but we have that, with positive probability, 
none of the edges in E2U E^U {e^} are updated from time 1 to time 2. Knowing this, 
again with positive probability, all edges in Ei are updated from time 1 to time 2 and all 
of them take exactly the value u := Zi(eo) (this is due to the Dirac mass du in the law 
Ue)- Since we started at equilibrium, Zt=2 has the equilibrium law, and edges in E2 are all 
open or all closed in the projections of Zt=2- This shows that with positive probability, at 
least N edges appear simultaneously as one raises p. 

Now, if A is getting very large, we can divide the box into a lattice of 2N x 2N squares. 
Starting from Za,o ~ /^A, the above strategy works in each box independently of what 
happens in other boxes. Stated like that, it looks wrong, since obviously the dynamics 
itself is not independent from one square to another, but all that is needed in the above 
procedure is a positive lower bound on the probability that this "scenario" happens. Using 
the structure of the dynamics, it is not hard to see that if yi, . . . , yx denote the indicator 
functions of the events that the scenario happened in the squares i £ { 1 , . . . , K} , then there 
is an independent product of Bernoulli e > variables which is stochastically dominated 
by our vector (yi, . . . ,yx). In particular, the emergence of clouds is somewhat ergodic in 
the plane. 

By changing slightly the argument, one can show that there are such clouds for any 
open interval of the variable p G [0, 1]. □ 

The above proof is not quantitative at all. Therefore, many natural questions on these 
emerging clouds remain: how do they look, how large are they, how does their law 
depend on the level p at which they appear? In particular, recovering the correlation 
length from such geometric considerations appears to be quite a challenging program. We 
shall discuss some of these questions in the next subsection. 



An intuitive explanation for the clouds. We end this subsection by a hand-waving 
argument why these clouds of simultaneously opening edges appear and may play an 
important role in the dynamics of any monotone coupling. Consider a monotone coupling 
(a;p, tJp+Ap). Due to the factor gi# clusters partition function, FK configurations Up 

tend to have as many clusters as possible. Without this factor, one would be in the case of 
q = 1, i.e., standard percolation, and the edge intensity would be exactly p. With q > 1, 
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the random-cluster configuration tries to maximize the number of clusters, hence the edge- 
intensity drops to a smaller value I{p) < p. In some sense, there is a fight between entropy 
(under the product measure edges ^-|^ — p)* edges ^ most configurations have edge- 

intensity p) and energy (which would correspond here to — log(g* clusters ^-j^ When one goes 
from p to p + Ap, new edges are added due to the entropy effect, but in such a way that 
not so many clusters will merge into a single one. A good strategy for adding many edges 
without a significant increase in energy is the following storing mechanism. Say we have 
two "neighboring" large clusters in Wp with closed edges going from one to the other (these 
closed edges are then large-scale pivotal edges). Once we decide to open one of them, it 
does not cost more energy to open a few others. 

Now, we have just seen that the monotone coupling is Markovian in p: in particular, 
the only way for this storing mechanism to actually happen is to have some values of p 
where the system can simultaneously open several edges. This indeed can happen, due to 
the atom in the update distribution, as shown in Lemma |3.8[ and the construction there 
was indeed a simple example of edges arriving simultaneously between two neighboring 
large clusters (the two components of £'2)- 

It is worth noticing that this heuristic explanation (based on entropy /energy consid- 
erations plus the Markov property) hints that this "non-linear phenomenon" should be 
much stronger near the critical point. Indeed, near Pc{q), there are many neighboring large 
clusters (i.e., many large scale pivotal points), which makes the storing mechanism more 
efficient. Away from criticality, this is not the case anymore. This intuition explains, for 
example, why we observe a blow-up of the derivative of the edge-intensity near pc, and 
why the emerging clouds are more important there. 



3.5 Questions on the structure of emerging clouds 

Finite volume case. The previous subsection shows that there are non-trivial clouds 
with positive probability. Let us consider the case of G = A„ := [— n, n]^ with free boundary 
conditions. (Another natural choice would be to consider discrete tori A„ := Z^/nZ^). We 
strongly suspect the following behavior: 

Question 1 (Macroscopic clouds near pc). For all n > 1, with ^\^-prohahility at least 
a universal constant c > 0, there is at least one macroscopic cloud in A^, i.e., whose 
diameter is larger than cn. Furthermore, with probability going to 1 as n ^ 00, the labels 
of such macroscopic clouds concentrate around the critical value pc{q = 2). 

To answer such a question, it is natural to run the dynamics at equilibrium (i.e., 
Zq ~ A*A„) for a short amount of time that is given precisely by the rescaling 

1 ^-13/24+0(1) _ 



n2^4(2) 



Doing so, only finitely many macroscopic pivotal edges will be resampled, and it is easy to 
convince ourselves that with positive probability at least two of them will pick the same 
label thus creating a macroscopic cloud. This intuition is close to being rigorous, since we 
have at our disposal a 'stability property' from the forthcoming |GP) which suggests that 
the "geometry" of ^tL^ could be recovered with high precision from Zq plus the updates of 
the initially macroscopically important edges (neglecting the "smaller" updates). However, 
the stability result holds only for ujp^, not the entire coupling Z. Thus, a certain control 
on the concentration of the labels around pc would be helpful for both parts of Question [TJ 
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The intuition that big clouds should appear only around the critical point can be 
translated into the following conjecture: 

Question 2 (Local clouds away from pc). For any 6 > 0, emerging clouds with labels 
outside of {pc — 5,pc + 5) are local in the sense that the largest such cloud in should he 
of logarithmic size. 

A natural way to attack this question would be via a couplin g arg ument. Namely, 
construct a coupling {Z^^"^^ , Zj^^"^^) (see the notation in Subsection 3.3) whose marginals 

are IJ^^~^^ ^ and whose coordinates are identical on a small neighborhood of the origin, but 
with probability at least A'^ (with A G (0,1)), are independent of each other outside a 
box of size k (an exponential decay of correlations). Such a statement is proved for the 
supercritical (or subcritical) random-cluster measure (pz^,p,qi Y^t the lack of a DLR (spatial 
Markov) property for our monotone coupling /iA„ prevents us for extending the result to 
the coupling in an obvious way. 

Finally, it would be interesting to prove quantitative results on the size of emerging 
clouds in the finite volume case (A„). This question is further discussed in |GH| . 



Infinite volume case. Clouds are well-defined objects: since there is a unique limiting 
measure 11^2 of Grimmett's monotone coupling and for any e S E{7?), cloud(e) can still 
be defined relatively to a sample Z from ^■^2. 

Question 3. Prove that a.s. there exist non-trivial emerging clouds. 

This does not follow directly from the existence of non-trivial clouds in the finite volume 
case. Assuming the above question, the next natural question would be the following: 

Question 4. Is it the case that emerging clouds are a.s. finite when q < 4:? 

Note that for large q, an infinite number of edges appear at Pc{q)- Again, |GH] contains 
some relevant discussion. 



3.6 What about the influence of an edge? 

To conclude this section, let us mention another natural approach to a "geometric" under- 
standing of the near-critical random-cluster model. 

As a continuation of the work by Kesten on near-critical percolation |Kes87j . Russo's 
formula should be replaced by a slightly different formula. Fix an increasing event A. As in 
the case of percolation, the intuition suggests that the derivative of <t^^Gpq^-^) with respect 
to p is mostly governed by the influence of one single edge, switching from closed to open. 
The following definition is therefore natural in this setting. The (conditional) influence on 
A of the edge e £ E, denoted by /^(e), is defined as 

^^(e) := (t>G,p,q('^\'^ is open) - (p^^^^^^^{A\e is closed). 
With this notation, we have the following formula: 

Proposition 3.9 (See |Gri06j ). Let q > 1 and e > 0; for any random- cluster measure 
4>G pq P ^ [^) 1 ~ ^^^y increasing event A, 

where the constants in x depend on q and e only. 
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It is tempting to use this extension of Russo's formula to see what our results on the 



correlation length (Theorem 1.2) may imply on the influences /^(e). To avoid boundary 
issues, let us consider the case of the torus T„ := Z^/nZ^, and let An be the event that 
there is an open circuit with non-trivial homotopy in T„. It is easy to check (by self- 
duality) that 0p^,2(^n) ^ 1/2- The results from Section |2] can easily be generalized to the 
torus. In particular, there exists a constant A > such that if pn ■= Pc(2) + A " , then 

</'p„,2(^n) >3/4. 



Pc 



V? 



Using the above Proposition 3.9 this says that 

(^l(e/^or) + /l(e^,er))dp>0(l)- 

where Chor and e^er are any horizontal and vertical edges in T„. Since it is natural to 
expect that on the interval [pcPc + A ^'°^" ], influences behave reasonably smoothly, the 
following conjecture should hold. 

Conjecture 3.10. For any n > I, X > 0, p £ [pc - \^^^^,Pc + A^^^p^] and any e £ T„, 

Il(e)>c- ^ 



n-y/logn 

where c = c(A) is some positive constant. 



In fact since it is reasonable to conjecture that in Theorem 1.2, one has actually 



I^P,e{p) ^ \p — Pc\ ^, one may strengthen the previous conjecture into the following one: 
Conjecture 3.11. For any n > 1, A > 0, p G [Pc — ^;Pc + ^] (ind any e G T„, 

c-</^ (e)<c-i-, 
n ^" n 

where c = c(A) is some positive constant. 

These conjectures are beyond reach with the techniques of the present paper. 

4 Other values of q 

As promised in the Introduction, most of this section will rely on predictions from physics 
to investigate what happens for general q G [1,4]. 

4.1 Some useful critical exponents 

Let us start by collecting several useful exponents. 

• ^1(9) denotes the one-arm exponent 

• £,4{q) denotes the four-arm exponent 

• a = Od{q) describes the behavior of the specific heat near Pc{<l)- That is, C^^''{p) ~ 
\p-Pc{qT^'''^. 
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P = /3(g) describes the behavior of the "magnetization": this can be interpreted as 
(j)p,q{0 ^oo)^{p- (as p \ pc{q)). 



v{q) corresponds to the correlation length: LF^''{p) ~ |p — Pc{(l)\ 



v{g)_ In 



• 7] = r](q) corresponds to the correlation function [x -f-T- y] ~ \x — y 

particular, assuming RSW, this exponent is twice the one-arm exponent (,i{q)- 

Let us summarize the physics predictions on these exponents in the following table. 
The expressions are simplified using the term u = u{q) = ^ arccos(^) = 2 — ^7^. 



Exponents 



u = u{q) 



a = a{q) 
V = v[q) 



6= 6(g) 



^4 = Uq) 



predictions 



arccos(^) 



2(1-2m) 
3(1-m) 

T+M 

12 
2—u 

3{l-u) 



2{2-u) 



4:{2-u) 
~3 



2-u 



Most of the previous critical exponents can be found, for example, in [Wu82j . except the 
four-arm or pivotal exponent which is more of a geometric nature. This latter exponent 
is computed in |DCG| using SLE^ calculations and assuming the (conjectured) correspon- 
dence 

An 

arccos(— ^) 

This SLE exponent was also derived by Wendelin Werner |Wer09a| . Assuming a proof of 
conformal invariance for the critical random-cluster model with parameter q as well as 
a proof of quasi-multiplicativity on the discrete level, this would say that 



FK, 



a. 



(n) 



n 



Such an estimate is of course far from reach at the moment, except for q = 1 (see |SW01) ) 
and q = 2 (see |DCGj ). but in this section we will assume that it holds. 

These exponents are not all independent of each other: there are scaling and hyperscal- 
ing relations between them that are expected to hold in large generality, across different 
models and underlying graphs. (But they are proved mathematically in very few cases 



only.) We will briefly discuss one of them in Subsection 4.3 



4.2 Near-critical behavior for q G (1,4] and self-organized monotone cou- 
pling 

Let us consider the random-cluster model on with fixed cluster-weight q £ (1,4] (we 
drop it for several notations). In this subsection, our goal is to illustrate that there is a 
strong self-organized mechanism within this monotone Markov process which goes beyond 
the specific heat effect. To show this, we will take for granted the specific heat exponent 
(which describes the blow-up of the derivative of the edge-intensity at criticality for q >2, 
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Figure 4.1: The blue curve corresponds to what the correlation length exponent in the 
critical random-cluster model with parameter q would be if the monotone coupling hap- 
pened to be "Poissonian". The red curve is a refinement of the blue one where the specific 
heat is taken into account. Finally, the black curve represents the actual correlation length 
exponent. 



as mentioned in Subsection 3.2), and based on this, we will estimate what would be the 
correlation length exponent if there was no self-organized mechanism (i.e., if new edges 
simply arrived in a Poissonian way). 

Let us first describe our setup: we will restrict ourselves to a finite but very large 
window An := [— n, n]^ with, say, wired boundary conditions. Now, starting from a critical 
random-cluster configuration ojp^ in A„, we raise p = pc + Ap until macroscopic effects start 
being non-negligible. If po is the value where we stop, we should thus obtain the relation 
L{po) X n. 

Now, the specific heat exponent a = a{q) has the following interpretation when q > 2: 



dp Mp-Pc(g)| 

where := 4>z'^,p,g{^ is open) is the random-cluster edge-intensity. From this result 

(which is at the level of a prediction in the physics literature), it is reasonable to expect 
that in the finite volume range, one has 

as far as n < L{p) (where In '(p) denotes the edge-intensity in A„ with, say, wired 
boundary conditions). This is known in the case q = 2 (where a{2) = with logarithmic 
blow-ups), where the finite-size behavior matches with the infinite volume one, as we 
already mentioned in Subsection |3.2[ 

On the other hand, when q G [^,2), the specific heat exponent a{q) G [—2/3,0) mea- 
sures "only" the second order variation of the partition function of the random-cluster model 
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around {p,q) = iPciq),Q), and there is no blow-up of the derivative of the edge-intensity 
in this case: 

dp 

as far as n < L{p). In particular, the specific heat exponent in our analysis will play a role 
only for q S [2,4]. (See, however, the next subsection.) 

Let us then start from a critical configuration LOp^ in A„ (with wired conditions) and 
let us raise p to the level p = Pc + in such a way that one still has n < L{p). From 
the above discussion, one expects that about Ap edges will arrive. If we 

assume the absence of self-organization, i.e., if edges arrive more or less independently of 
the current configuration (except possibly a local rate which would depend on whether 
the endpoints of the edge are connected or not), then each of these arrivals should be 
macroscopic pivotal flips with probability about 

n-MQ) (we implicitly harnessed the fact 
that the pivotal exponent does not vary below the critical length). Therefore, at n ~ L{p), 
we expect 

Let LP™""(p) denote the correlation length obtained via the above analysis. We find 



which is represented as a function of q by the red curve in Figure 4.1 (the blue curve 
represents the result of the same analysis when specific heat blow-ups are not taken into 
account). 



As one can see from Figure 4.1 the actual correlation length L(p) is much smaller than 
^Poiss^^^ when q G (1,4], which reveals for all these random-cluster models non-trivial 
self-organized schemes as p increases near Pc{q)- 

4.3 The hyperscaling relation between correlation length and specific 
heat 

There is the following well-known hyperscaling relation between correlation length and 
specific heat: 

2- a = ud, (4.1) 

which is expected to hold for all q, provided that the dimension d of the underlying lattice 
Z"^ is low enough. 

We have been arguing that the mere quantity of new edges arriving does not explain 
the correlation length alone, but the self-organized structure in which they arrive also 
matters. Nevertheless, the hyperscaling relation tells us that the correlation length is in 
fact determined by the change in the edge intensity, in some other way. Although we have 
not managed to relate the mechanism for the hyperscaling to self-organization, let us give 



a brief heuristic derivation of (4.1 ) in the case of g > 2 (where it is linked to the change in 
the edge intensity), in a way easily accessible to mathematicians. 

On the infinite lattice, the derivative of the edge intensity is dl{p) /dp k, \p—pc\~°^- Since 
a finite system starts looking similar to the infinite system at the scale of the correlation 
length L{p), it is reasonable to expect (as we did also in the previous subsection) that the 
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edge intensity is changing in a similar way in a finite system of side-lengtli L{p). Hence, 
tlie number of new edges in this box, when getting from pc to p, should be about 



1— Q |„ „ \ud+l—a 



^[L{p)]d,pH-^[L{p)]d,pA^\ = Hp) \p-Pc\ " ~ \p-p 



coming from integrating the above rate of change. On the other hand, similarly to (3.8) 



we have ^^[L(p)]d^p \^\ ^ Varj^j-p^jd ^ which means that the standard deviation of |a;| 
from its mean should be around \p—pc\^~'^'^~"')^'^. The near-critical window should be given 
by the near-critical bias just starting to overcome the existing fluctuations, which results 
in the equation 

—i/d — a 
—ud + l — a = , 



yielding the desired identity (4.1). 



For the percolation case {q = 1), Section 9.2 of |Gri99| contains a similarly vague 
heuristics, while |BCKS99] studies those hyperscaling relations that do not involve the 
specific heat. 



5 On the dynamical sensitivity of random-cluster models 

We end this paper with a slightly tangential discussion on the heat-bath dynamics for the 
random-cluster models with q > 1. A rigorous treatment of the special case q = 2 will be 
presented in [GP]. 

The purpose of the first subsection is to highlight the following interesting phenomenon: 

Phenomenon 5.1 (conjectural). There exist "natural" critical two-dimensional systems 
(i.e., scale-invariant and so on) with the property that they have pivotal points at all scales 
(hence are expected to be noise sensitive), but for which there are no exceptional times of 
infinite clusters along the natural heat-bath dynamics (hence are not dynamically sensitive). 

The second subsection is a rigorous treatment of the case q > 25.72 (which should then 
also hold for all g > 4 but is yet conjectural). 



5.1 Noise- versus dynamical sensitivity for q E [1,4] 



The (conjectural) Phenomenon 5.1 is rather surprising, since, in every model examined 
in the literature so far, noise sensitivity and dynamical sensitivity {i.e., the existence of 
exceptional times with an infinite cluster) have been present or absent together. One can 
see this as another illustration of the fact that a good understanding of the structure of 
the set of pivotal points is by far not enough to answer questions on noise- or dynamical 
sensitivity. Recall that all the current proofs of noise- and dynamical sensitivity for critical 
percolation (g = 1) rely on the "largeness" of the Fourier spectrum of percolation, a very 
different object compared to the pivotal points (see [GPSlOa] ). The Fourier spectrum of 
random-cluster models should still be "large", i.e., concentrated on "high frequencies" once 
projected on an orthogonal basis of eigenfunctions of the heat-bath dynamics. Nevertheless, 
it seems that, when q gets closer to 4, the Fourier spectrum is not large enough to overcome 
the typical critical behavior (the smallness of the one-arm probability), hence exceptional 
times will not exist. We now discuss the reasons for this phenomenon in greater detail. 

As for the dynamics for the monotone coupling, an infinite-volume heat-bath dynamics 
on can be constructed for every {p,q) with q > 1. These dynamics are constructed 
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Figure 5.1: The red curve represents the upper bound on the Hausdorff dimension of 
exceptional times for random-cluster models, q G [1,4] which is obtained heuristically in 
this subsection. It highhghts an interesting transition as the variable q increases. Note 
that the red curve gives the correct bound of 31/36 for q = 1 (proved in |GPS10a] ). 

as monotone limits of the wired or free finite volume heat-bath dynamics and are unique 
as soon as the infinite-volume measure is unique at {p, q) (even though, theoretically, 
they could differ from each other). Fix g S [1, 4] and p = Pc{q)- In this case, the infinite- 
volume measure is expected to be unique. Set (wj '')t>o to be the unique critical dynamics 
obtained either from the free or the wired limit. The natural question raised by the 
discussion of the percolation case is the following: 

Question 5. For q £ (1,4], are there exceptional times t almost surely for which there is 
an infinite cluster in ' ? // "yes" and if £q denotes the random set of these exceptional 
times, what is the (almost sure) Hausdorff dimension of £q? 

We conjecture the following property for which we will then give a heuristic proof based 
on the above exponents from the physics literature. 

Conjecture 5.2. Let 



q* :=4cos^(-V14) w 3.83. 

For q £ {q*,4], FKg percolation is NOT dynamically sensitive. 
If q £ [l,q*], one has a.s. 



- Suiqr - 8n(g) ' ^'''^ 



where, as above, u{q) := ^ arccos(^) = 2 — - 



(9)- 
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This conjecture hints that there is a critical qc G above which there are no 

exceptional times. It is natural to expect that q i— t- E [dim((?|j)] is continuous and thus that 



> 1. It is quite possible that qc = q* and (5.1) is an equality. 



Heuristic explanation of the conjecture. For standard percolation, the usual way 
to prove upper bounds on the dimension of exceptional times (and prove sometimes their 
non-existence) is to divide the time evolution into tiny intervals / of length e each, domi- 
nate the union of critical configurations {cjf : t G /} by a slightly supercritical percolation 
configuration, and use the near-critical behavior. However, due to the complicated struc- 
ture of the near-critical FKg percolation highlighted in the previous section, one cannot 
easily do this for q > 1. Instead, one can rely on a more artificial near-critical version, 
where new edges are added to a critical random-cluster configuration in a Poissonian way. 
The correlation length L{p) of this near-critical model can be derived from the four- arm 
critical exponent as in the case q = 1 (this assumes RSW and the analog of near-critical 



stability described in Phenomenon 1.5 and Subsection 3.1). This gives 

1 



L{p) 



\P-Pc{q)\ 



2-«4{'j) 



Below this correlation length, the artificial near-critical configuration has the same con- 
nectivity properties as a critical random-cluster configuration. In particular, if all our 
configurations uj^ , t G I are dominated by a configuration u}p^+e, one obtains 

p[/n£:, /0] <p[o^^Z(p, + e)] 

<F[O^L{pc + e)] 

Since one needs 0(e^^) intervals / to cover [0,1], a first moment argument implies that 
a.s. 

'""''<^'> ^ ' - ^ ■ 

Plug ging in the expressions from Subsection |4.1| explains the second part of Conjec- 



ture 



solution 



5.2 



For the first part, it is enough to solve the equation 2^£^'l{q) = 1) so that above its 
', the dimension would be "negative", i.e., we expect the set £q to be almost surely 



empty. The solution of this equation is given by q* . 



Finally, let us mention that (5.1 ) will be made rigorous in the special case g = 2 in the 
forthcoming paper |GP) . 

5.2 Random-cluster models with g > 4 

The purpose of this subsection is to briefly explain what occurs when g > 4. In this 
case, the phase transition is expected to be of first order (proved for q large enough, i.e., 
q > 25.72): in ]p^^«'^'''"^ there is no infinite cluster, while in jp^^'J^^i^'ed |-j^gj.g Another way 

^ ' Pc{q) ' Pc(q) ^ 

to formulate this phenomenon is that, in a finite system [— n, n]*^, the critical parameter 
Pcisi) of the infinite system is already outside the critical window: it is subcritical in the 
free case, supercritical in the wired case. However, the finite system is not more interesting 
even at the "true n-dependent critical value" (with any reasonable definition; say, the left- 
right crossing probability is exactly 1/2): any infinite- volume measure should be a mixture 
of and Pp (again, proved in any dimension d > 2 for large enough q, 
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see |LMMS^9l j and the references there). In particular, the correlation length remains 
uniformly bounded for all n G N and p G [0, 1]. Therefore, most of the previous discussion 
of the near-critical regime does not make sense. 

What one could still ask is the size of the critical window, and how it is related to the 
structure of the emerging clouds in the monotone coupling. It appears that the clouds are 
even larger than in the case q G (1,4]. First of all. Theorem 4.63 and Proposition 8.59 of 
|Gri06] say that there is an atom at p in the infinite volume monotone coupling if and only 
if there is non-uniqueness of the measure at that p. That is, in the infinite volume limit, a 
positive fraction of the edges appear at Pc{q) for q large. (One can also prove this from the 



exponential decay of connectivity in the critical free measure, see (5.2) below.) A priori, 
for the finite systems this means only that there is an interval shrinking to Pc in which a 
positive fraction of the edges appear. However, the uniform boundedness of the correlation 
length strongly suggests that the transition from the subcritical to the supercritical phase 
can happen only by a simultaneous appearance of a positive density of edges; i.e., there 
should be a uniformly positive atom even for the finite systems. 

One can also discuss the possibility of noise and dynamical sensitivity. Yet, the expo- 
nential decay of correlations suggests that the model is not noise sensitive. As an illustra- 
tion, we exploit the known exponential decay estimates for 1?^^^)^'^'^'' for q > 25.72 in order 
to obtain the following result. Of course, the theorem is expected to hold for all g > 4. 

Theorem 5.3. When q is large enough, there are no exceptional times for the infinite free 
boundary heat-bath dynamics on critical FKg configurations. 

Proof. The proof is based on |Gri06| Theorem 6.35] which states that in dimension d = 2, 
and if q > 25.72, then at the critical point Pc{q) = Psdid), one has for the free infinite 
volume limit: 

^M,)A^^d[-n,n]') < C exp(-c(g)n) , (5.2) 

where c{q) > is a positive constant which depends only on q > 25.72. By monotonicity, 
this result implies that for the random-cluster measure on the finite box A^r := [—NjN]"^ 
endowed with free boundary conditions, then for all radius n such that 2n < N, one has 

<^L,Pc(g),g(^^^(^ + ^")) ^ C exp{-c{q)n) , 

for all points x G A„. This in turn implies that for all x G A„, the probability to have 
a four-arm event around x of radius n is bounded above by the same exponential bound 
exp(— c(g)n). 

Now, let us fix a large radius n ^ 1. Our goal is to find a small upper bound on 

for the free boundary infinite- volume heat-bath dynamics on Z^. Since, as we discussed 
earlier, this infinite volume limit is obtained as an increasing limit of finite volume heat- 
bath dynamics and since the event under consideration is a cylinder event, one has 

g{n) = ^lim^,/.o^^p^(^)^^(3t G [0, 1] s.t. ^ dAn) . 

Using the random variable Xn = X^^^^ to denote the number of flips over the time 
interval [0, 1] for the event {0 i — > 9A„}, one easily obtains 

E[X„] < 0{l)n? sup </'i^,p,(g),,(x ^ X + d{K/2)) 

x€An 

< 0{l)n^ exp(-c(g)n/2) . 
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This in turn implies the bound 



<f>LMi),i(^* ^ [0' 1] ^ ^^") ^ Oil) r? exp(-c(g)n/2) , 

which gives a uniform upper bound in N. In particular, g{n) < 0(1) exp(— c(g)n/2) 
and thus, taking n — >■ oo, one concludes that a.s. there are no exceptional times for the 
critical free random-cluster measure on 7? with q large enough. □ 
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